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In  these  lecturesAseveral  methods  are  presented  that  are  useful  for 
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PREFACE 


The  David  W.  Taylor  Lectures  were  initiated  as  a  living  memorial  to 
our  founder,  in  recognition  of  his  many  contributions  to  naval  architecture 
and  naval  hydrodynamics.  Admiral  Taylor  was  a  pioneer  in  the  use  of  hydro- 
dynamic  theory  and  mathematics  for  the  solution  of  naval  problems.  The 
system  of  mathematical  lines  developed  by  Taylor  was  used  to  develop  many 
ship  designs  for  the  Navy  long  before  the  computer  was  invented.  He 
founded  and  directed  the  Experimental  Model  Basin;  perhaps  most  Important 
of  all,  he  established  a  tradition  of  applied  scientific  research  at  the 
"Model  Basin"  which  has  been  carefully  nurtured  through  the  decades  and 
which  we  treasure  and  protect  today.  In  the  spirit  of  this  tradition,  we 
invite  an  eminent  scientist  in  a  field  closely  related  to  the  Center's  work 
to  spend  a  few  weeks  with  us,  to  consult  with  and  advise  our  working  staff, 
and  to  give  a  series  of  lectures  on  subjects  of  current  interest. 

Our  third  lecturer  in  this  series  is  Professor  John  V.  Wehausen  of 
the  University  of  California,  Berkeley.  Professor  Wehausen  gradurted  from 
the  University  of  Michigan  and  received  his  Ph.D.  in  mathematics  from  there 
in  1938.  He  is  no  stranger  to  the  Center  or  its  problems.  From  1946  to 
1949  he  worked  here  as  a  mathematician  and  developed  his  interest  in  naval 
hydrodynamics.  He  later  became  the  resident  mathematician  in  the  Depart¬ 
ment  of  Naval  Architecture  at  the  University  of  California,  where  he 
further  developed  his  knowledge  of  the  theory  of  water  waves  and  the  hydro¬ 
dynamics  of  ships.  Professor  Wehausen  has  trained  many  outstanding  re¬ 
searchers  in  the  field  of  naval  hydrodynamics  through  his  teaching  at  the 
University,  and  we  are  most  honored  that  he  agreed  to  be  a  David  W.  Taylor 
lecturer. 
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ABSTRACT 

In  these  lectures  several  methods  ere  presented  that 
are  useful  for  solving  free-surface  boundary-value  problems: 
separation  of  variables  and  the  Fourier  method,  the  method 
of  reduction  and  reflection,  the  method  of  Green  functions, 
and  the  method  of  multipoles.  Each  method  Is  illustrated 
by  one  or  more  examples.  In  the  examples  the  fluid  Is  assumed 
to  be  inviscid  and  incompessslble  and  the  flow  irrotational. 
The  boundary  conditions  have  been  linearised.  The  examples 
themselves  are  all  concerned  with  diffraction  and  forced 
motion.  Although  the  methods  are  applicable  to  a  much  wider 
class  of  problems,  this  restriction  allows  a  simple  formula¬ 
tion  of  the  physical  problem  and  immediate  involvement  with 
the  method  itself. 


INTRODUCTION 

The  analysis  of  an  engineering  problem  in  fluid  mechanics  usually 
proceeds  along  the  following  lines.  First  one  selects  a  model  of  a  fluid: 
Navier-Stokes  or  inviscid,  compressible  or  not,  with  or  without  surface 
tension,  etc.  The  decision  can  usually  be  based  upon  an  examination  of 
various  physical  parameters  characterizing  the  problem.  Next  a  de¬ 
cision  must  be  made  with  regard  to  the  flow:  is  it  steady  or  unsteady, 
is  it  irrotational,  etc.?  Again,  this  will  be  based  upon  physical 
parameters  (or  perhaps  a  desperate  need  to  avoid  complications).  The  next 
step  is  an  exact  formulation  of  the  equations  of  motion  and  of  the  bound¬ 
ary  and  initial  conditions  for  the  fluids  and  bodies  present.  (Ideal¬ 
ly  this  should  be  accompanied  by  existence  and  uniqueness  studies.) 

The  problem  will  usually  be  too  difficult  mathematically,  so  the  next 
step  is  to  replace  these  equations  by  other  simpler  ones.  The  region  o.' 
usefulness  of  these  simpler  equations  will  also  depend  on  certain 
parameters  associated  with  the  problem.  Finally,  assuming  that  the 
equations  are  now  tractable,  one  proceeds  to  solve  them  either  analyti¬ 
cally  or  perhaps  by  direct  numerical  approximation  methods. 


One  must  grant  Immediately  that  the  distinction  between  these  steps  is 
not  always  clear*  Is  approximating  the  Navier -Stokes  equations  by  the  Euler 
equations  really  different  from  approximating  the  flow  about  an  elongated 
body  by  that  given  by  a  slender-body  approximation?  One  might  argue  that 
the  first  involves  a  physical  parameter  associated  with  the  fluid  and  the 
latter  a  geometrical  parameter,  but  both  are  approximations.  As  further 
cases  are  considered,  one  realises  that  the  boundaries  are  fusty. 

However,  these  problems,  important  and  fascinating  as  they  are,  will 
not  concern  us  here  for  we  are  going  to  limit  ourselves  to  just  one  aspect 
of  the  procedure  outlined  above,  namely,  the  finding  of  analytical  solutions 
once  a  tractable  problem  has  been  formulated.  Indeed,  we  shall  restrict 
our  attention  still  further,  namely,  to  problems  with  a  free  surface  in 
which  the  fluid  is  assumed  to  be  inviacld,  incompressible,  and  subject  to  a 
gravitational  force  and  in  which  the  flow  is  assumed  to  be  lrrotational. 

In  addition  to  these  restrictions  we  shall  also  assume  that  the  boundary 
conditions  are  linearised. 

It  is  obvious  that  we  are  bypassing  a  host  of  interesting  questions. 

We  shall  not  consider  how  or  whether  the  linearised  equations  we  deal  with 
can  be  made  part  of  a  systematic  approximation  scheme.  This  problem  will 
be  relegated  to  one  of  the  other  "steps"  described  above.  We  have  re¬ 
stricted  the  purview  of  these  lectures  so  drasticcally  for  the  following 
reasons.  Although  a  single  topic  cuuld  be  selected  and  developed  fairly 
Intensively,  in  a  set  of  eight  lectures  each  lecture  would  have  to  be 
built  upon  the  preceding  ones.  The  audience  would  necessarily  have  to 
devote  some  thought  between  lectures  to  the  topic  under  discussion  in  order 
to  keep  the  whole  development  in  mind.  And,  of  course,  regular  attendance 
would  be  necessary.  Since  the  situation  is  different  from  a  university, 
where  both  outside  study  and  regular  attendance  can  be  expected,  it  seemed 
more  useful  to  select  a  subject  that  can  easily  be  divided  into  one-  or 
two-lecture  units  so  that  not  too  much  baggage  has  to  be  accumulatd  as  one 
proceeds.  I  believe  and  hope  that  this  will  not  make  the  lectures  less 
interesting. 
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The  lectures  will  consist  of  s  discussion  of  the  severel  methods  In 
use  for  finding  snelytlcal  solutions  for  free-surfece  problems*  Each 
method  will  be  illustrated  by  one  or  more  examples ,  which  ve  hope  will 
themselves  have  some  intrinsic  interest*  Several  of  the  examples  could 
be  solved  by  more  than  one  of  the  methods,  and  it  would  undoubtedly  be 
Interesting  to  do  this.  The  original  intention  was  not  to  do  so  in  order 
to  obtain  a  greater  diversity  of  physical  problems.  However*  the  examples 
treated  do  tend  to  fall  into  one  area:  diffraction  and  forced  harmonic 
motion.  This  was  not  necessary.  It  just  so  happened  that  such  problems 
seemed  to  be  among  the  simplest  with  which  to  illustrate  the  various 
methods. 

The  methods  will  be  chiefly  the  following:  separation  of  variables 
and  the  Fourier  method,  the  method  of  reduction  and  reflection,  the 
method  of  Green  functions,  the  method  of  multipole  expansions,  and,  time 
permitting,*  variational  methods.  I  shall  try  to  Indicate  some  of  the 
advantages  and  disadvantages  of  each  and  the  limitations  in  their  use. 

Since  the  easiest  problems  have  usually  been  solved  first,  the  examples 
will  often  not  be  from  the  recent  literature.  On  the  other  hand,  they 
are  also  the  most  useful  for  purposes  of  illustration  just  because  they 
are  fairly  simple.  I  shall  try  to  compensate  for  this  by  calling  attention 
to  recent  literature  in  which  the  method  under  discussion  has  been  used. 

I  do  not  wish  to  imply  that  solving  more  complex  problems  simply  required 
the  investigators  to  turn  the  crank  a  little  harder.  Difficulties  almost 
always  arise  in  applying  one  of  the  methods  to  a  new  problem,  and  it  may  be 
necessary  to  alter  it  in  some  appropriate  way. 


It  did  not. 


PART  1:  THE  BASIC  EQUATIONS 

As  coordinate  system  we  begin  with  an  inertial  one  fixed  in  the  fluid. 
We  choose  it  rectangular  and  right-handed  and  with  Cy  directed  oppositely 
to  the  force  of  gravity  in  both  two  and  three  dimensions.  The  undisturbed 
water  surface  will  be  taken  as  the  (x,z)  plane.  As  stated  above,  we  suppose 
the  fluid  to  be  inviscid  and  incompressible  and  the  motion  irrotational. 
Consequently,  there  exists  a  velocity  potential  $  (x,y,z,t)  with  absolute 
velocity 

V  =  grad  $  =  V$  =  (1.1) 

Conservation  of  mass  gives 

A$  =  4>  +$  +  $  “  0  (1.2) 

xx  yy  zz 

The  three  conservation-of-momentum  equations  (Euler  equations)  reduce  to 
Euler  integral 

•j  O 

$t+  2  |V*|  +  gy  +  p/p  =  constant  (1.3) 

If  the  free  surface  is  represented  by 

y  =  Y (x,  z,  t) 

then  the  kinematic  boundary  condition  on  the  free  surface  is 

4>  (x,Y,z,t)Y  -  *  +$>Y  +  Y  =  0  (1.4) 

x  x  y  zz  t 

and  the  dynamic  boundary  condition  is 

*t(x,Y,z,t)  +  \  |V$|2  +  gY  =  0  (1.5) 
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If  a  solid  boundary  is  present,  then  the  kinematic  boundary  condition  states 
that 


$  =  v*n 

n 


(1.6) 


at  each  point  of  the  surface,  where  V  is  its  velocity  at  that  point.  If 
the  surface  is  given  by 


F(x,y,z,t)  =  0 


this  condition  may  be  written  in  the  form 


$  F  +  $  F  +4>F  +  F  =  0 
xx  y  y  z  z  t 


(1.7) 


(1.8) 


Dynamic  boundary  conditions  may  also  be  necessary  on  solid  boundaries. 

For  example,  in  the  case  of  freely  floating  bodies,  these  take  the  form  of 
the  Newton  equations  for  motion  of  a  rigid  body. 

We  have  stated  earlier  that  we  shall  consider  only  linearized  problems 
without  discussing  the  rationale  of  the  linearization.  The  Laplace  equa¬ 
tion  is  already  linear,  but  the  solution  of  a  nonlinear  problem  will  not, 
of  course,  be  the  same  as  that  of  a  corresponding  linearized  problem.  The 
two  free-surface  conditions  become 


<fy(x,0,z,t)  -  Yt  =  0 


$t(x,0,z,t)  +  gY  = 


If  Y  is  eliminated,  this  condition  becomes 


(1.9) 

(1.10) 


$>  (x,0,z,t)  +  g$>  =  0 

tu  y 


d.ii) 


Linearization  of  the  boundary  conditions  on  solid  boundaries  is  more  com- 
licated  for  it  will  depend  upon  the  nature  of  the  problem  and  the  basis 
of  the  linearization. 


.  »  j 
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The  free-surface  condition  (1.11)  may  appear  In  different  forms, 
depending  upon  the  circumstances.  If  the  motion  is  assumed  to  be  harmonic 
in  time  with  frequency  0,  we  may  write 

$>(x,y,z,t)  -  ^(x.y.z)  cos  at  +  $2(x,y,z)  sin  at 

-  Re  [<t>e~i0t],  ♦  -  +  i<{>2  (1.12) 

Then  (1.11)  becomes 

4>(x,0,z)  -  ^  <j>  -  0  (1.13) 

a  y 

If  it  is  convenient  to  use  a  moving  coordinate  system  x^x+ct,  y=y, 
z  ■  z,  and  if  <|>(x,y,z,t)  =  $(x,y,z,t)  then  (1.11)  becomes 

c2(f> _ (x,0,z,t)  -  2c<j)_  +  <|>tt  -  c<j>_  +  g(j>_  *  0 

xx  xt  x  y 

If  the  motion  is  steady  in  the  moving  coordinate  system,  then  <j>  does  not 
depend  on  t  and  the  three  middle  terms  drop  out. 

We  assume  that  these  equations  and  their  provenance  is,  infact,  known. 
They  are  given  only  for  ready  reference. 


PART  2:  THE  FOURIER  METHOD 

Thi3  is  probably  the  oldest  and  most  widely  used  method  for  solving 
boundary -value  problems  in  mathematical  physics.  It  is  the  classical  method 
of  separation  of  variables,  determination  of  eigenfunctions,  and  expansion 
of  the  solution  in  them.  The  expression  may  take  the  form  of  either  a  sura 
or  an  integral,  and  we  shall  give  an  example  of  each. 

A  limitation  of  the  method  in  application  results  from  the  fact  that 
we  usually  would  like  to  identify  the  surface  upon  which  boundary  values 
are  imposed  with  coordinate  surfaces  of  the  variables.  This  explains  why 
so  many  of  the  examples  of  classical  physics  deal  with  rectangular,  circu¬ 
lar,  and  elliptical  cylinders  and  with  rectangular  solids,  spheres,  ellip¬ 
soids,  etc.  These  all  fit  into  some  system  of  variables  in  which  the 
equations  are  separable.  However,  the  possibilities  are  limited  for  there 
are  only  eleven  coordinate  systems  in  which  the  three-dimensional  Laplace 
equation  can  be  separated  and  five  in  which  the  two-dimensional  one  can  be 
separated  (see  e.g.,  an  article  by  N.  Levinson,  B.  Bogert,  and  R.M.  Redheffet 
on  pages  241-262  of  Volume  7  (1949)  of  the  Quarterly  of  Applied  Mathematics 
and  references  cited  there) .  We  shall  be  concerned  only  with  coordinate 
systems  in  which  one  of  the  coordinates  is  y. 

SEPARATION  OF  VARIABLES 

We  shall  look  for  a  solution  in  the  form 

$(x,y,z,t)  =  X(*.«)  *(y)  T(t)  (2.1) 


Substitution  into  A$  *  0  yields 


(Xxx  +  XM>  YT  +  XY"T 


0 


or,  if  wa  assume  x  and  Y  are  not  identically  zero. 


v< 


X  Y 


(2 


* 


7 


X 


T 


Each  term  on  the  left-hand  side  must  be  a  constant,  say,  a  and  $,  respec¬ 
tively,  where 

a  4 J  *•  0  (2.3) 

Hence 


A2X  -  ax  -  0,  Y"  -  BY  -  0  (2.4) 

It  was  stated  above  that  one  of  the  situations  in  which  application  of 
separation  of  variables  is  most  convenient  .13  where  one  of  the  boundaries 
is  a  coordinate  surface.  The  surface  y  ■  0  is  already  such  a  surface.  Let 
us  choose  the  bottom  to  be  another  such.  We  shall  consider  two  cases.  In 
one,  a  horizontal  bottom  exists  at  y  ■  -h.  In  the  other,  the  fluid  is 
infinitely  deep.  In  the  first  case,  the  boundary  condition  (1.6)  becomes 

<Jy(x,-h,z,t)  »  0  (2.5a) 

In  the  second  case,  we  shall  relax  this  to  the  condition 

|$y(x,y,z,t) |  <M<“asy+-“  (2.5b) 

For  the  function  form  (2.1),  these  imply 

Y'(-h)  °  0  and  |  Y  *  (y )  J  <N<“asy->-“  (2.6) 

We  next  substitute  (2.1)  into  the  free-surface  condition  (1.11)  to 
obtain 


X(x,z)l(0)T"(t)  +  gxY’T  -  0 


T"(t)  +  8  Txo)  T  ■  0 


We  shall  restrict  attention  to  the  case  where  the  coefficient  of  T  Is  posi¬ 
tive  and  write 


a2  -  gY'(0)/Y(0) 


(2.7) 


The  solution  for  T  is  then 

T  *  A  cos  at  +  B  sin  at  (2.8) 

Let  us  now  consider  the  equation  for  Y, 

Y"  -  0Y  =  0 

2 

If  we  choose  0  -  k  >0,  then 

Y  =  C  eky  +  D  e"ky 

The  boundary  conditions  in  (2.6)  then  yield 

Y  ■  cosh  k(y  +  h)  and  Y  *  eky  (2.9) 

respectively.  (The  multiplicative  constant  can  be  safely  dropped  since 
this  is  already  provided  for  in  T  in  (2.8).) 

But  now  (2.7)  gives 

2  2 

a  *  gk  and  O  ■  gk  tanh  kh  (2.10) 


for  h  *  00  and  h  <  00 ,  respectively. 


From  (2.3)  and  (2.4),  we  then  find  that  x  must  satisfy 


a2X  +  fc2X  ■  0  (2.11) 

2 

Next  we  choose  3  ■  -k  <0.  Then 

Y  *  C  cos  ky  +  D  sin  ky 

The  boundary  condition  Y'(-h)  =  0  yields 

Y  -  cos  k(y  +  h)  (2.12) 

For  infinite  depth,  boundedness  imposes  no  further  condition.  However, 
from  (2.7)  we  deduce 


a 


2 


so  that  we  may  write 

2 

Y  =  cos  ky  +  —  sin  ky  (2.13) 

Equation  (2.13)  reflects  only  the  effect  of  the  free  surface  and  not  of  the 
bottom. 

Equation  (2.12)  substituted  into  (2.7)  gives 


<?  -  -  gk  tan  kh 


(2.14) 


The  equation  for  x  that  is  associated  with  3  *  -  k^  is 


V>X  "  k  X  “  0 


(2.15) 
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Let  us  now  suppose  that  a  is  fixed  and  ask  what  values  of  k  will  be 
determined  for  h  <  ®  by  (2.10)  and  (2.14).  This  is  easily  seen  from  a 
graphical  display.  We  write  (2.10)  and  (2.14)  in  the  form 

2  2 

a  h  1  .  ,  , ,  .  oh  1 

7  kh*  tanh  kh  ™d  7“  th  “  "  tan  kh  (2.16) 

A  graphical  solution  of  the  first  equation  shows  a  single  solution  (see 
Figure  1). 


kh 


Figure  1  -  Graphic  Solution  of  First  Equation  in  (2.16) 

By  the  same  method  one  finds  an  infinite  number  of  solutions  for  the 
second  equation  (see  Figure  2) . 

For  the  case  h  <  ®,  we  have  now  obtained  the  following  solutions  for  #: 

XQ(x,z)  cosh  kQ(y  +  h)  [AQ  cos  ot  +  BQ  sin  at] 

2  (2.17) 

a  -  gkQ  tanh  kQh,  A2x  +  kQ2x  -  0 


ll  k) 
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and 


X^x,*)  cos  k^y  +  h)  [A^  cos  at  +  sin  at3 


2  2 
a  ■  -  gk^  tan  k^h,  A2X^  ”  "  0 


(2.17) 


A  sum  of  such  terms  (appropriately  convergent)  will  evidently  also  be  a 
solution.  With  regard  to  the  family  of  functions 

{cosh  kQ(y  +  h),  cos  k^(y  +  h),  cos  k2(y  +  h) . . . . }  (2.18) 


it  is  possible  to  prove  the  following  property:  Any  pair  of  functions  is 
orthogonal  on  the  interval  (~h,0),  i.e., 


(2.19) 


J  cosh  k«(y  +  h)  os  k . (y  +  h)  dy  ■  0 
-h 


t  ~ k‘ 


(y  +  h)  cos  k.  (y  +  h)  dy  ■  0  if  i  s'  j 


One  nay  prove  also  the  following  theorem: 

Theorem:  If  f(y)  is  defined  and  square- integrable  on  the  interval  [-h,0), 
chen 


f(y)  r  bQ  cosh  kQ(y  +  h)  +  bj,  cos  k^y 

i-1 


+  h)  (2.20) 


where 


brt  - 


0  sinh  2k 


4k  r 

ST+20  J.  ‘'V>  00811  Vv  +  h>  d> 


O  ■'O  -h 


and 


4  k, 


sin  2kjh  +  2k ^h 


f°  * 

-h 


(y)  cos  k^y  +  h)  dy 


The  formulas  for  the  b's  follow  from  the  orthogonality. 
For  h  ■  we  have  the  following  solutions  for  $: 
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(2.21) 


X(x,x)  eR/  [A  cos  ot  +  B  ain  ot], 
a2  -  gk,  A2x  +  k2x  -  0 

X(x,*)  J^cos  ky  +  ^-  gin  ky"|  [A  coa  ot  +  B  sin  ot], 
k  arbitrary,  A2X  *  k‘‘x  ■  0 

Do  we  have  anything  analogous  to  the  representation  theorem  (2.20)Y  There 
is,  in  fact,  an  analog,  first  given  by  Havelock  in  1929,  that  is  reducible 
to  the  Fourier- integral  theorem.  It  is  as  follows: 

Theorem:  If  f(y)  is  absolutely  integrable  on  (-  0),  then 


,®  ,0 


i  {  dk  l 


dn  f(n)  C-k  cos  ky  +  v  sin  ky)  (k  ccs  kn  +  V  sin  kn) 

k2  +  v2 


(2.22) 


+  2v  e' 


J  f(n)  evn  dn,  v  -  o2/g 


We  may  write  this  in  a  form  more  closely  analogous  to  (2.20)  by  pre¬ 
senting  it  as  a  pair  of  transforms: 


f (y)  *  bg  e  y  +  J  b(k)  (k  cos  ky  +  v  sin  ky)  dk 


bQ  -  2V  J  f (n)  e  dn 

•■OQ 


(2,23) 


u 

1  f 

TT  J 


k  cos 


2  2 
k  +  v 
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f(n)  dn  co®  k(y  -  n),  with  f(-n) 


f(n) 


-  f(y> 


by  the  usual  Fourier  integral  theorem.  We  have  made  use  above  of  the  fol¬ 
lowing  two  integrals: 


cos  mx 

2  2 
a  +  x 


dx 


TT 

2a 


x  sin  mx 


2  . 

a  i*  x 


2 


,  TT  - 

dx  ■  j  e 


m  >  0,  a  >  0 


Further  separation  of  variables  should  now  depend  upon  the  nature  of 
the  problem  to  be  solved.  In  one  that  we  shall  consider  below,  a  wavemaker 
at  one  end  of  a  rectangular  channel,  it  is  natural  to  use  again  rectangular 
coordinates  and  to  assume  x(x»*)  “  X(x)Z(s).  If  one  wished  to  find  the 
wave  motion  generated  by  an  oscillating  vertical  post  of  circular  cross 
section,  one  would  use  polar  coordinates.  However,  as  mentioned  earlier, 
the  number  of  possible  configurations  is  very  limited. 

We  shall  determine  the  functions  X(x)  and  Z(s)  for  separation  in  rec¬ 
tangular  coordinates.  Substitution  into  (2.11)  quickly  gives 


(2.24) 


or 


X"  -  ctX  -  0,  z"  -  yZ  -  0,  a  +  y  ■  <  0 


If  y  -  q  >0,  then 


Z  -  C  eq*  +  D  e~q* 


and 


(2.25) 


X  -  E  cos  (k2  +  q2)1/2x  +  F  sin  (k2  +  q2)1/2x 
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If  Y  “  -  q  <  0,  than 


Z  -  C  cos  q*  +  D  sin  qs  (2.26] 

and 

X  -  E  cos  (k2  -  q2)1/2x  +  F  sin  (k2  -  q2)l/2x  if  q2  <  k2 

X  -  E  axp  (q2  -  k2)1/2x  +  F  exp  £-(q2  -  k2)1/2xj  if  q2  >  k2 

If  we  substitute  into  (2.15),  we  again  obtain  (2.24)  except  that  now 
ot  +  y  ■  k2.  Now  if  y  ”  q2  >  0* 

Z  «  C  eqi  +  D  e‘qE  (2.27] 

and 

X  -  E  cos  (q2  -  k2)1/2x  +  F  sin  (q2  -  k2)1/2x  if  q2  >  k2 

X  -  E  exp  (k2  -  q2)1/2x  +  F  exp  £-(k2  -  q2)1/2x]  if  q2  <  k2 

2 

If  Y  “  ~q  <  0,  then 


Z  »  C  cos  q*  +  D  sin  qs  (2.26! 

and 

X  -  E  exp  (k2  +  q2)1/2x  +  F  exp  £-(k2  +  q?)1/2xj 

This  is  a  complete  census  of  the  possibilitites  in  rectangular  coordinates 
Let  us  now  turn  to  some  applications.  The  first  one  is  a  wavemaker 
at  one  end  of  a  semi-infinite  rectangular  channel.  The  second  is  a  wave- 
maker  in  a  wall  bounding  the  region  x  >  0,  y  <  0.  The  treatments  parallel 
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each  other  except  that  the  y-  and  z-dimensions  are  bounded  in  one,  and  in¬ 
finite  in  the  other.  The  first  employs  Fourier  series  in  the  solution, 
the  second  Fourier  integrals. 

WAVEMAKER  IN  A  CHANNEL 

Consider  a  semi- inf inite  channel  bounded  by  plane  walls  at  y  ■  -h  and 
z  ■  0  and  z  ■  b.  At  x  *  0  there  is  situated  a  wavemaker  that  moves  accord¬ 
ing  to  the  equation 

x  *  F (y, z)  sin  at  (2.29) 

The  channel  extends  to  infinity  in  the  direction  Ox.  In  practice,  this 
can  be  approximated  by  having  a  very  efficient  wave  absorber  at  the  er.d 
away  from  the  wavemaker.  A  wide  variety  of  wavemakers,  both  two-  and  three- 
dimensional,  can  be  described  by  a  proper  choice  of  F(y,z).  We  3hall  fur¬ 
ther  assume  that  the  motion  has  persisted  for  a  long  time,  so  that  transient 
motions  associated  with  starting  the  wavemaker  have  died  out  and  the  fluid 
motion  is  also  harmonic  with  frequency  a. 

Let  us  formulate  the  mathematical  problem.  We  may  evidently  take  $  in 
the  form  (1.12)  and  the  free-surface  condition  in  the  form  (1.13).  The 
boundary  conditions  on  the  sides  and  bottom  will  be 

4>  (x.y.O)  =  <j>  (x,y,b)  -  0,  q>  (x,  -h,  z)  -  0  (2.30) 

«  «  y 

The  linearized  boundary  condition  on  the  wavemaker  corresponding  to  (1.6)  is 

4>x(0,y,z,t)  =  aF(y,z)  cos  at 
or 

<t>lx(°»y»z)  *  aF»  *2x(0’y’z)  "  0  (2.31) 

There  is  still  one  missing  condition,  but  it  will  be  more  interesting  to 
see  this  forced  upon  us  later  on. 
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In  order  to  solve  this  problem,  we  shall  try  to  use  the  family  of 
solutions  given  in  (2.17)  together  with  the  solutions  for  x  given  in 
(2.25)-(2.28).  First  of  all  we  note  that  the  boundary  conditions  on  the 
side  walls  given  in  (2.30)  cannot  both  be  satisfied  with  any  combination 
of  the  exponential  solutions  (2.25)  or  (2.27).  These  conditions  will  be 
satisfied  by  the  trigonometric  solutions  (2.26)  and  (2.28)  if  we  take  C 
and  D  such  that 


Z  ■  cos  z,  m  ■  0,1,2,.... 

D 


(2.32) 


In  choosing  the  X's  to  go  with  the  Z's,  we  must  remember  that  in  (2.26) 
k  must  be  k^  and  in  (2.28)  k  may  be  any  one  of  the  k^  (see  (2.17)).  Finally, 
since  we  are  interested  only  in  bounded  solutions,  we  must  discard  the  in¬ 
creasing  exponential  in  the.  exponential  solutions  for  X.  If  we  now  sum  a]  1 
these  elementary  solutions  with  arbitrary  multipliers,  we  anticipate  that 
the  solution  to  the  problem  can  be  written  iu  the  following  form: 


4>(x,y,z) 


v2al/2 


.  .2 

11/2 

(b  ) 

J  * 

WTT 


cosh  kg(y  +  h)  cos  z 


b„  [-  ((r  )  -  k02  )  *]  oosh  k0(y  +  h) 


miT 

cos  i^—  z 


m=M+l 


mir 


cos  k^(y  +  h)  cos  -g—  z  (2.33) 


mir 


where  M  is  the  largest  integer  m  such  that  r- -  <  kn 


This  function  satisfies  the  Laplace  equation,  the  free-surface  condition, 
and  the  conditions  on  the  bottom  and  sides  of  the  channel;  4*  may  be  either 


or  4>2* 

Let  us  examine  the  condition  (2.31)  on  the  wavemaker: 

M  2  l/2 

4>x(°»y»z)  m  bm  [k02  ’(?)  ]  C08h  k0(y  +  h)  cos  “  Z 

m*0 

00  2  1/2 

+  JT  ~  bm  [( if)  “  kQ2  ]  cosh  k0(y  +  b)  cos  f-  z 
m-M+1 

00  00  _  2  1 1/2 

+  £ bim  [ki2  +  (?)  ]  COS  ki(y  +  h)  COS  ?  Z  (2*3A) 
1“1  m=0 


If  we  now  use  the  orthogonality  of  the  functions  {cosh  kQ(y  +  h) , 

cod  k^(y  +  h) ,  i  =■  1,2,...}  Stated  in  (2.19)  and  the  well-known  orthogonality 

of  the  family  {cos  z,m  **  0,1,2,...}  on  the  interval  0  <  z  <  b,  we  may 

D 

easily  derive  the  following  formulas  for  the  b's: 


—  2  1/2 

8^  (Binh  2k„h  +  2k„h)  [k02  -  (f  )  ]  bm 


U  D 

J  dy  J  dz  4'x(0,y,z)  cosh  kQ(y  +  h) 


-h  “0 


b 

8kr 


nrr  .  „ 

cos  t—  z,  m  <  M 

D 


(slnh  2k0h  +  2k0h)  [(f)'  -  k,,2]^2  bm 


(2.35) 


J  *J 

- h  0 


dz  4>x(C,y,z)  cosh  kg(y  v  U)  cos  ^  z,  m  >_  M  +  1 
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(sin  2^  +  2^) 


1/2 


b 


m 


dz  4>x(0,y,z)  cos  k^y  +  h)  cos  -g-  z. 


i  >  1,  n  >  0 


In  order  to  determine  <f>^,  we  replace  4>x  in  (2.35)  by  aF;  to  determine 
</>2»  by  0.  As  long  as  is  such  that  it  is  never  equal  to  mv/b  for  some 
integer  m,  all  coefficients  b  are  uniquely  determined.  The  completeness 
of  the  two  orthogonal  families  for  representing  functions  on  -h  £  y  £  0 
and  0  £  z  _<  b  tells  us  that  any  square-integrable  function  F  can  be  so 
represented.  For  <j>2>  the  b's  are  obviously  all  zero. 

Have  we  now  determined  the  functions  <j>^  and  Evidently  not,  for 

a^,  ....  aM  are  not  determined  for  either  4^  or  <j>2>  Since  we  have  now  satis¬ 
fied  all  our  boundary  conditions,  we  must  conclude  that  either  the  problem 
does  not  have  a  unique  solution  or  else  that  we  have  not  completely  formu¬ 
lated  it.  Here  the  latter  is  true.  If  there  had  been  another  wall  at 
x  »  i  >  0,  we  should  have  had  another  boundary  condition  to  determine  the 
coefficients.  We  evidently  need  something  to  replace  it.  This  is  the 
radiation  condition,  which  states  that  waves  must  propagate  down  the  chan¬ 
nel,  i.e.,  in  the  direction  Ox.  The  terms  exponential  in  x  play  no  role  in 
this  condition  for  they  die  out  as  x  increases,  representing  only  a  local 
disturbance  near  the  wavemaker.  In  order  to  apply  the  radiation  condition, 
we  may  write  the  solution  in  the  complete  form 

$(x,y,z,t)  ■  <t>1  cos  at  +  4>2  sin  at 

and  then  choose  the  a's  so  that  each  of  these  first  M  terms  represents  a 

progressive  wave  moving  to  the  right.  It  is  easy  to  see  that  this  is 

achieved  by  taking  a  ^  -  0  and  a  ^  *  b  where  the  superscripts 

m  mm 
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A, 

/ 


'  ) 


t 


respectively  refer  to  and  c? , ,  for  then  the  first  M  terms  take  the  form 


M 


ZV*  -”{h2-(f)] 


2  t  1/2 


X  - 


m^O 


mir 


at  }  cosh  k^(y  +  h)  cos  z 


This  choice  of  the  a's  has  been  forced  upon  us  by  the  radiation  condition, 
the  latter  being  necessary  to  achieve  a  unique  solution.  However,  if  the 
problem  had  been  formulated  as  an  initial-value  problem  in  which  the  motion 
started  from  rest,  this  solution  would  have  been  obtained  automatically. 

It  will  be  convenient  for  discussion  of  its  properties  to  write  the 
solution  in  a  slightly  different  form.  Define 


8k 


b  [sinh  2k^h 


0  f  f 

h  '-T 2k~hf  J  dyJ  dz  F(y>z)  COsh  k0(y  +  h)  COB  IT  Z 


■oh]  -h  "0 


(2.36) 


8k 


0 


"im 


b  [sin  2k^h 


i+~2k~hT  Jh  dyjf  dz  cos  kj^y  +  h>  cos  2 


Then  $  is  given  by 


M 


<i>(x,y,z,t)  ■27 


[W-(r)T* H 


-  [*o2-(r)] 


2-il/2  m 


C  sin 


mTT 


x  cosh  k^(y  +  h)  cos  -g—  z 
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inril'ilki  iJi  ;  ii  ii.— I,--,.  :  k.Ul <t’llWi!fr¥Tif 


o 


x  cos  t—  z  cos  Ot 
b 


00  00 


yy _ 2 _ 

1-1  w-0  [k±2  +  (fl)  ] 


,2-,  1/2 


27172  Cim 


-  [kj2  +  x  cos  k± (y  +  h) 


miT 

x  cos  -g—  z  cos  Ot 


(2.37) 


According  to  (1.10),  the  free  surface  itself  is  given  by  Y  ■  -g-1<I>t(x,0,z,t) 
It  is  as  follows: 


Y(x,z,t)  -  kQ  sinh  kQh  (r  2  .2-.  1/2  \ 

-[•.■-m1]1'  *[  r  1  ' 

-z  y-1..  i-[(F),-.-]u,.)-f- 

m-M+l  ^  -  k  2  1  f 


z  sin  Ot 


00  00 


+  y  y ——-r\n  c,  -p{-  k2  +  (?)Y/24 

w  To k2 -((ai) ]  i  Ll  u;j  I 


mir  . 

cos  z  sin  ot 

D 


(2.38) 


Let  us  examine  the  solution.  It  is  evident  that  it  breaks  down  for 
any  kg  such  that  kg  *  mir/b  for  some  integer  m.  Indeed,  the  formulas  deter¬ 
mining  the  bQ  in  (2.35)  do  not  do  so  for  the  particular  m  for  which  this 
happens,  so  that  we  have  not,  in  fact,  found  a  solution  for  such  a  value 


ip>  a»'»n  -  ■■ 


of  kg  (or  its  associated  frequency  a).  Hence,  there  is  an  infinite  sequence 
°l*a2’°3* *  *  *  ^or  w^ich  no  solution  appears  to  exist.  Let  us  examine  the 
behavior  of  (2.38)  in  the  neighborhood  of  these  frequencies. 

We  shall  make  a  thought  experiment  in  which  the  frequency  a  starts 
out  quite  small  and  then  increases  little  by  little  through  these  special 
values  (Figure  3). 


Figure  3  -  Critical  Frequencies  for  Wavemaker  in  a  Channel 

When  a  is  quite  small,  kgti  <  nh/b  and  M  ■  0.  After  the  local  dis¬ 
turbance  has  died  out,  there  is  then  a  single  two-dimensional  progressive 
wave  propagating  down  the  channel: 

Cq  sinh  kgh  cos  (kgX  -  ot) 

However,  as  a  approaches  from  below,  k^h  -*  irh/b  and  the  coefficient  of 
the  first  exponential  term  in  the  second  summation 
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grows  unboundedly.  At  the  same  time  the  coefficient  of  x  In  the  exponen¬ 
tial  d-  reases,  so  that  the  disturbance  no  longer  dies  out  so  quickly. 

If  o  is  Increased  just  beyond  o^,  then  M  ■  1  and  the  first  summation 
consists  of  two  terms 


and 


Cq  sinh  k^h  cos  (kgX  -  at) 


;g  sinh  kgh 


'-(?)] 


2  -i  1/2 


cos 


(k-a)  ] 


2  t  1/2 


x  -  at 


IT 

COS  z 


The  first  is,  if  course,  a  two-dimensional  progressive  wave.  The  second 
represents  a  progressive  wave  of  longer  length  sloshing  from  side  to  side 
as  it  progagai.es.  It  might  be  indicated  schematically  as  shown  in  Figure  4. 
If  a  is  close  to  a^,  its 


• 

CREST  : 

• 

• 

• 

• 

e 

• 

•  • 

•  • 

;  TROUGH  : 

•  • 

a  • 

• 

Figure  4  -  Representation  of  Sloshing  Progressive  Wave 


amplitude  will  be  very  large.  As  a  increases,  its  amplitude  will  decrease 
and  approach  sinh  kgh.  However,  before  0  Increases  very  much,  it  will 
approach  and  another  "catastrophe”  will  occur.  After  passing  Og,  a 
third  tvu -  of  progressive  wave  modulated  by  cos  (2ir/b)z  will  be  added  to 
the  first  uwo  (s^  Figure  5). 
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Figure  5  -  Representation  of  Third  Type  of  Progressive  Wave 


It  is  evident  that  every  time  0  passes  through  a  value  Om,  there  is 
an  associated  crisis  during  which  a  new  progressive  wave  is  added,  begin¬ 
ning  with  a  very  large  amplitude.  These  frequencies  are  sometimes  called 
"cutoff  frequencies."  One  should  note  that  even  though  the  wavemaker  is 
carrying  out  a  three-dimensional  motion,  the  resulting  progressive  wave 
will  be  two-dimensional  if  a  is  small  enough,  i.e. ,  a  <  o^.  On  the  other 
hand,  if  the  motion  of  the  wavemaker  is  exactly  two-dimensional,  so  that 
F(y,  z)  ■  F(y),  then  Cm  »  0  for  m  ^  1  and  also  oim  ■  0  for  m  ^  1.  Thus, 
none  of  these  crises  occurs.  However,  if  the  supposed  two-dimensional 
wavemaker  is  only  a  little  out  of  true,  one  may  anticipate  the  catastrophic 
behavior  described  above. 

One  may  argue  legitimately,  that  the  linearized  theory  Is  no  longer 
a  valid  approximation  in  the  neighborhood  of  these  cutoff  frequencies. 
However,  the  behavior  described  above  does  occur  and  is  well  known.  It  is 
not  the  same  as  a  similar  phenomenon  known  as  "cross-waves,"  recently 
studied  by  Garrett  (1970) ,  Mahony  (1972) ,  and  others  (Garrett  gives  a  his¬ 
tory  of  this  phenomenon) .  Cross-waves  are  a  nonlinear  phenomenon.  I  am 
not  sure  who  first  solved  the  problem  we  have  just  examined.  It  is  included 
in  the  Biesel  and  Suquet's  (1951)  encyclopedic  article  on  wavemakers,  and 
certain  aspects  have  been  studied  in  more  detail  by  Kravtchenko  (1954) . 
However,  it  is  accessible  to  anyone  familiar  with  the  methods  of  classical 
mathematical  physics. 
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One  question  may  remain.  Since  we  have  not  found  a  solution  for  the 
cutoff  frequencies  o^,  what  should  be  done  about  them?  These  are  like  the 
resonance  frequencies  in  forced  motion  of  a  harmonic  oscillator.  One  must 
reformulate  the  problem  as  an  initial-value  problem.  Then  one  may  antici¬ 
pate  that  the  solution  for  these  frequencies  will  grow  without  bound  as 
t  -*>. 


WAVEMAKER  IN  A  WALL 

As  an  example  of  the  use  of  Fourier  integrals  and  of  the  representa¬ 
tion  theorem  (2,22),  (2.23),  we  have  chosen  a  problem  similar  to  the  one 
just  discussed.  The  chief  difference  is  that  the  bottom  and  side  walls 
are  removed  so  that  the  only  boundary  is  the  wall  at  x  *  0,  The  fluid  is 
in  the  region  x  >  0,  y  <  0.  The  wall  Itself  is  flexible  and  moves  accord¬ 
ing  to  the  equation 


x  ■  F(y,z)  sin  Ot 


(2.39) 


where  we  assume  that 


|F(y,«)| 


<  oo 


(2.40) 


As  before,  we  assume  that  the  motion  has  become  harmonic  in  time. 

The  mathematical  formulation  is  similar  to  that  of  the  wavemaker  in  a 
channel.  We  take  ♦  in  the  form  (1.12)  and  the  free-surface  condition  in 
the  form  (1.13).  The  first  two  conditions  of  (2.30)  are  abandoned  and  the 
third  replaced  by 

lim  <t>  (x,y,z)  -  0  (2.41) 

V  —  03  ' 
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The  boundary  condition  (2.31)  still  hoLds.  In  addition,  here,  as  before, 
we  need  an  explicit  condition  stating  that  the  waves  propagate  away  from 
the  wavemaker. 

We  must  now  make  a  selection  of  the  solutions  (2.21),  (2. 25)- (2. 28) . 
The  exponential  solutions  in  z  must  be  discarded  because  they  become  un¬ 
bounded  on  one  side  or  the  other.  The  decreasing  exponentials  in  x  are 
allowable.  The  most  general  solution  for  $  that  satisfies  these  conditions 
is  the  following: 


-v 

<l>(x,y,s)  -J  dq  | 
0  { 


2  2  1/2 

cos  [(V  -  q  )  x]  [AQ(q)  cos  qz  +  BQ(q)  sin  qz] 


-  sin  [(v2  -  q2)^2x  ]  [CQ(q)  cos  qz  +  DQ(q)  sin  qz]  j 

■*' /  dq  exp  [-(q2  -  V2)^2x]  (CQ(q)  cos  qz  +  DQ(q)  sin  qz]  e^ 


QO  CO 

+  J  dkj  dq  exp  [-(q2  +  V2)^2x]  [C(k,q)  cos  qz  +  D(k,q)  sin  qz] 
0  0 


x  [k  cos  kz  +  v  sin  ky] 


(2.42) 


Here,  <p  can  be  either  or  *  We  can  confirm  immediately  that  (2.41)  is 

satisfied  if  the  integrals  exist.  For  the  single  integrals,  this  follows 

Vy 

from  the  form  of  the  factor  e  .  For  the  double  integral,  this  follows 
from  the  Riemann-Lebesque  Lemma. 

Let  us  now  compute 
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00 

.-.u 


.2  .2,1/2 


b0<«>  *  “JQ  dq  “  q4|**  lC0<q>  coa  qt  *  D0(q^  8in  qs^ 


.  -0 


■  2v  J  ^(O.y.x)  dy 


and 


b(k,s)  -  - J  dq(q2  +  k2)1^2  [C(k,q)  cos  qx  +  D(k,q)  sin  qx] 


r° 

2  |  ^  v  k  cos  ky  +  \>  sin  ky  . 

■  ?  J  — rr— j — 1  d” 


k  +  v 


The  ordinary  Fourier  inversion  integral  now  gives 


0  00 


- |„2  -  ,v/2 19,(0  ♦  1D0(,)1  -  it  J  *  J  *•  ♦,»•»•*>  •vy  '+1,‘ 

—  (2.45) 


and 


0  00 


-  (k2  +  q2)1/2  [C(k,q)  +  iD(k,q)l  -  J  dy  J  dx  *x(0,y,x) 


k  cos  ky  +  v  sin  ky  iqx 
k2  +  v2  * 


The  coefficients  A^(q)  and  Bg(q)  have  not  yet  been  determined  for 

either  ^  or  ^  For  $2*  since  ^2x^,y'*^  "  a**  the  C  s  and  D'8  are 

xero  and  $2  takes  the  form 
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*2(x»y**>  "  J  **  co*  Hv2  -  q2)1/2x)  cos  q.  + 


02  cos  q t  +  b02  *in  q«l  « 


vy 


If  v  consider  only  the  part  of  ^  cos  ot  +  ♦  .in  ot  involving  the  fir.t 
integral,  we  find 


J  dq  (cos  I(v2  -  q2)1/2x]  {AQl  cos  q«  + 


Bqi  sin  qa]  coa  ot 

2  2  1/2 

sin  [(v  -  q  )  x]  [Cq^  cos  qc  +■  Dq^  sin  q«]  coa  ot 
+  cos  t(v2  -  'i2)1/2»l  t*02  cos  „  +  b02  sin  ,t)  sin  ot  J  s'* 
In  order  that  this  should  represent  an  outgoing  wave,  we  must  set 


A01  "  B01  ?  °*  A02  “  +  C01*  b02  *  +  D- 


01 


This  gives  for  this  integral 


'[  “1  8l“  «V2  -  lV/2*  -  ot)  [C„,  co.  o,  +  n 


01  cos  qz  +  D01  Sin  qzj 


2  |  dq  {C01  sin  l(v2  "  +  q*  -  ot J  +  CQ1  sin  [(v2  -  q2)1/2} 

«qa  -  Ot)  -  Dq1  cos  t(v2  -  q2)1/2x  +  qz  -  ot)  +  DQ  sin  [(v2  -  q2)1^2* 


-  q*  -  Ot)  )e 


) 


yy 


(2.46) 
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v  / 


The  functions  C^,  Dq^,  C,  and  D  are,  of  course,  known  from  (2.45)  with 
4>x  replaced  by  aF(y,z). 

The  solution  has  no  especially  interesting  properties  like  those  occur¬ 
ring  in  the  channel.  This  fact  is  associated  not  with  the  infinite  depth 
but  with  the  infinite  width.  The  cutoff  frequencies  would  also  have  occurred 
here  if  there  had  been  walls  at  z  »  0  and  z  ■  b.  In  this  case  the  solution 
would  have  been  a  bit  different  in  that  we  would  have  used  Fourier  series  in 
the  z-direction  rather  than  the  Fourier  integral. 

Let  us  compute  the  average  rate  at  which  work  is  being  done  by  the 
wavemaker.  According  to  the  linear  approximation,  the  instantaneous  rate 
is 


W(t) 


0 

I  *>  \ 


m  0  00 

dz  p (0,y , z)  u(0,y,z)  =  ~p  J  dyj  dz  $>x(0,y,z,t)  ^t(0,y,z,t) 
—00  —00  —00 


■  -  pjTdydz  [<j;^x  cos  at  +  <j>2x  sin  at]  a  [-<^  sin  at  +  (j>2  cos  at] 

■  -  po  JJdydz  [<l,lx<(,2  co®2  at  ~  ^2x^1  8in2  0t 

+  ^1x^1  +  ^2x^2^  sln  Dt  COS  at^ 

■  -  pa2  JJ dydz  F(y,z)  [<j>2  cos2  at  -  <j>^  sin  at  cos  at] 

where  we  have  made  use  of  the  boundary  conditions  for  <J>^x  and  <{)2x  when  x  •  0. 
If  we  wish  to  find  only  the  average  rate,  we  must  still  integrate  with  re¬ 
spect  to  t  over  one  period  2tr/a.  This  yields 


1 


W  -  PO2  JJ  dydz  F(y,z)  MO.y^) 


-  +  -pa  2  JJ  dydz  F(y,z)  J 


dq  tCQ1  cos  qz  +  Dni  sin  qz]  e 


vy 


01 


-  +  |pO  2  JJ  dydz  F(y, z)  J  dq  [CQ1  +  iDQ1]  e"iqz  eVy 
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Using  (2.45)  with  <p^x(Q,y,z)  «  0F(ysz),  we  find 


v 

W  -  -J-  P03J]  <iy<lz  F(y,z)  Jo  dq  ^ 


evy  e'1<lz 


^JJdndq  F(n,C)  evr| 


0  00  0  00  -V 

ru  P  i~  J  dy  I  dz|  dn  I  dS  F(y,z)  F(n,C)  J  dq  ~ eV(y+n) 


—00  —00  —00  —00 


0  00 


(V  -  q  ) 

r2  P  f~  I  dy  J  dz  |  dn  j  F^y*z>  F<n,C)  eV(y+n)  JQ  [v(z  -  C)1 

(2.47) 


—CO  —00  —00  — oo 


This  formula  is  somewhat  analogous  to  the  Michell  integral  for  the  resis¬ 
tance  of  a  thin  ship.  Indeed,  the  latter  problem  is  a  kind  of  steady-state 
analog  in  which  one  considers  the  flow  past  a  small  (s  "thin")  bump  in  a 
wall.  In  fact,  Michell  derived  this  wave-resistance  formula  by  a  Fourier- 
integral  method  very  similar  to  that  just  used;  it  is  described  in  Wehausen 
(1973,  pp.  143-148). 

The  problem  discussed  above  was  first  solved  by  Havelock  (1929).  As 
he  pointed  out,  (2.22)  or  (2.23)  can  be  used  together  with  any  other  coor¬ 
dinate  system  in  which  one  can  separate  the  function  x(x»z)  (2.21).  He 
applies  the  method  to  the  waves  generated  by  the  oscillation  of  a  vertical 
circular  post,  but  he  considers  only  the  simplest  case  when  there  is  no 
dependence  upon  the  angles.  An  oscillation  with  angular  dependence  would 
lead  to  situations  of  the  sort  encountered  with  the  wavemaker  in  a  channel. 

The  problem  we  have  considered  is  essentially  a  simple  one.  However, 
the  method  can  be  and  has  been  applied  to  more  complicated  situations.  For 
example,  Ursell  (1947)  has  used  it  to  Investigate  diffraction  of  water  waves 
from  an  obstacle  in  the  form  of  a  vertical  plate  of  finite  length,  treating 
it  as  a  two-dimensional  problem  with  normal  incidence  on  the  plate.  Evans 
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and  Morris  (1972)  have  extended  the  treatment  to  oblique  incidence*  Such 
problems  require  finding  a  solution  on  each  side  of  the  plate  and  then 
matching  them  at  the  interface  below  the  plate.  In  addition,  one  must  fake 
into  account  conditions  at  the  edge  of  the  plate.  A  method  that  can  be 
used  to  solve  a  diffraction  problem  can  also  be  used  to  solve  a  related 
forced-motion  problem. 


PART  3:  THE  REDUCTION  METHOD 

The  essence  of  this  method  is  to  consider  at  first  the  function 

instead  of  ♦  itself.  This  function  is  also  a  solution  of  the 

tt  y 

Laplace  equation  and  vanishes  on  the  surface  y  «  0.  Let  us  define 

H(x,y,z,t)  -  <j»tt  +  g<j)y  (3.1) 

Ue  extend  the  region  of  definition  of  H  from  those  parts  of  the  lower  half- 
space  y  <  0  occupied  by  fluid  to  their  mirror  Images  in  y  >  0  by  the  equation 

H(x,-y,z,t)  *  -H(x,y,z,t)  (3.2) 

Since  H(x,-0,z,t)  «  0,  it  follows  immediately  that  not  only  is  H  continuous 

on  y  *  0  but  also  H  ,  H  ,  H  and  H  .  However,  it  also  follows  from  the 

definition  that  Hy(x,-0,z,t)  *  Hy(x,+0, z,t)  so  that  H^  is  also  continuous. 

Furthermore,  since  H  *»  -H  -H  ,  it  follows  that  also  H  is  continuous. 

*  yy  xx  zz’  yy 

Hence,  H  has  been  extended  as  a  harmonic  function  into  the  upper  half- 
spacn.  If  the  rest  of  the  boundary  conditions  can  be  carried  over  to  H  and 
do  not  complicate  the  problem,  the  solution  for  H  may  turn  out  to  be  easier 
than  that  for  <&.  Of  course,  after  finding  H,  we  are  still  confronted  with 
the  problem  of  finding  4>,  which  may  not  be  easy. 

Although  the  method  is  presented  above  in  a  three-dimensional  formula¬ 
tion,  it  has  found  its  greatest  use  in  two-dimensional  problems  where  the 
powerful  methods  of  analytic-function  theory  may  be  applied.  We  make  a 
small  digression  to  review  some  of  the  basic  facts  concerning  irrotational, 
two-dimensional  flow  of  an  incompressible  fluid. 

The  equations  expressing  conservation  of  mass  and  irrotationality 

u+v-0,  u-v-0  (3.3) 

x  y  y  x 

are  just  the  Cauchy-Riemann  equations  for  u  and  -v,  and  consequently 
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w  «  u  -iV  is  an  analytic  function  of  z  *  x  +  iy.  Furthermore,  the  first 

equation  ir.  (3.3)  implies  that  there  exists  a  stream  function  Y(x,y,t)  with 

¥  ■  -V,  'V  ■  u.  But  then  t  and  V  also  satisfy  the  Cauchy-Riematm  equations 

x  y 

and 

F(z,t)  -  4>(x,y,t)  +  mx.y.t)  (3.4) 

is  also  an  analytic  function.  It  is  easy  to  see  that 

F'(z,t)-w  (3.5) 

We  shall  call  F  the  complex  potential  and  w  the  complex  velocity. 

Let  us  now  consider  the  combination 

H  (z)  -  Ftt  +  igF'  -  *tt  +  i4'tt  +  ig(Ox  -  i*y) 

-  *tt  +  g*y  +  iWtt  +  gty)  (3.6) 

It  follows  from  the  free-surface  condition  that 

i  • 

Re  H(x  +  10)  -  0  (3.7) 

But  then  H(z)  can  be  extended  as  an  analytic  function  from  the  domain 
y  <  0  to  the  domain  y  >_  0  by  the  Schwarz  reflection  principle: 

H(x  -  iy)  -  -  H(x  +  iy)  (3.8) 

Tula  of  course,  is  the  analogue  of  (3.2)  and  as  in  the  three-dimensional 
case,  we  hope  to  be  able  to  exploit  this  extension  to  the  whole  plane  to 
find  an  easier  solution  for  H. 
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The  method  has  been  used  by  so  many  persons  that  I  shall  not  try  to 
sketch  its  history  except  to  say  that  1  believe  It  was  first  used  in  water- 
wave  problems  by  Levi-Civita.  It  has  been  extensively  exploited  by  various 
Russian  hydrodynamicists  but  also  by  others  in  the  United  States  and  else¬ 
where.  To  illustrate  the  method,  we  shall  use  diffraction  from  a  vertical 
plate  immersed  into  the  water  to  a  depth  l.  This  problem  was  treated  in 
this  way  by  Haskind  (1948) .  It  can  be  treated  equally  well  by  the  method 
of  Fourier  integrals,  and  in  fact  this  was  done  by  Ursell  (1947).  The 
method  has  recently  been  applied  by  Evans  (1970)  to  the  problem  of  diffrac¬ 
tion  about  a  completely  submerged  flat  plate.  It  is  almost  obvious  that  if 
the  diffraction  problem  can  be  solved  by  this  method,  then  the  problem  of 
the  waves  generated  by  small  oscillations  of  the  plate  can  also  be  solved. 
This  is  also  Included  in  the  cited  paper  by  Evans.  Later  on  we  shall  come 
back  to  the  reduction  method  when  we  treat  the  method  of  multipole 
expansions. 

DIFFRACTION  OF  WAVES  ON  A  VERTICAL  PLATE 

We  suppose  that  a  flat  plate  is  immersed  to  a  depth  l  in  an  infinitely 
deep  fluid  and  is  subjected  to  oncoming  waves.  Some  part  of  these  will  be 
reflected  and  some  part  transmitted,  as  indicated  schematically  in  Figure  6. 
In  the  neighborhood  of  the  plate  there  will  be  a  local  disturbance  that  we 
have  not  tried  to  represent. 
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The  incident  wave  will  be  represented  by 


Yj  •  A  cos  (kx  +  Ot) ,  a  ■  gk 


(3.! 


with  velocity  potential 


♦T  -  -  eky  sin  (kx  +  ot) 
I  o 


-  ^  eky  sin  kx  cos  at  -  ■§—  eky  cos  kx  sin  at 
a  o 


<PIC  cos  at  +  <J>IS  sin  at 


(3.H 


The  associated  complex  potential  of  (3.4)  is  easily  verified  to  be 


Fr(z)  «  -i  e_ikz  cos  at  -  ^  e“ikz  sin  at  (3.1 


-  fIC  cos  at  +  fIS  sin  at 

We  shall  denote  the  diffracted  wave  by  $D  so  that  the  velocity  poten 
tial  for  the  total  motion  itt 

♦  <t>  cos  at  +  <t>  sin  at 

I  D  Tc  s 

-  (<t>lc  +  $DC)  COS  ot  +  (<{>IS  +  4>ds)  Sin  at  (3.1 

Analogously,  the  complex  potential  is 

F  *»  F_  +  F_*  f  _  cos  at  +  fc  sin  Ot 

1 
1 

] 

-  (f__  +  f^J  cos  ot  +  (f  g  +  fns)  sin  at  (3.13)  j 


The  following  boundary  conditions  oust  be  satisfied  on  the  flat  plate: 
*x(0.y»t)  ■  0,  or  <j>cx(0,y)  -  $sx^0*y*  "  °*  0  -  y 

On  the  free  surface,  (1.11)  must  be  satisfied  by  4  and  4^  and  hence  (1.13) 
by  4C,  4S>  4  ,  and  4dS*  In  complex  form  these  are  the  following: 

On  the  plate: 

R  f’(0  +  iy)  -  Re  f '  (0  +  iy)  -  0,  0  >  y  >-2.  (3.15) 

6  C  S 

or  in  terms  of  FQ, 

Re  {  \o  eky  +  f^,(0  +  iy)  J  -  Re  f^.(0  +  iy)  -  0  (3.16) 

On  the  free  surface: 


Im  {f *  +  ikf (x  +  10)}  -  0 


where  f  is  any  one  of  f  .  f  ,  fDC,  fDg. 


(3.17) 


In  addition,  4Q  must  represent  outgoing  waves  at  a  distance  from  the 
plate.  In  complex  notation  this  may  be  expressed  as  follows: 


lim 

x  -*•  +  « 


(fDC  -  kfDSJ 


0,  It.  (t;s  +  kf*.) 

X  *►  +  ® 


(3.18) 


It  is  also  necessary  to  specify  the  behavior  of  the  flow  near  the 
sharp  edge  at  (0,  -  4)  for  there  will  be  a  singularity  at  this  point.  In 
order  to  limit  the  power  of  this  singularity,  we  shall  suppose  that 
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| z  +  Uj1/2  |f  *  <*)  I  <  M  near  z  -  -U 


or  in  ocher  words  that 


(z  4-  U)1/2  f’(z)  (3.19) 

Is  regular  at  Che  point  -  IS..  This  restriction  Is  In  conformity  with 
experience  at  similar  cusps  in  hydrodynamics. 

Finally  we  impose  some  further  conditions  at  infinity: 

| f  *  (z)  |  <  B  if  h|  >  l  +  1,  y  <  0  (3.20) 


and 


lim  f'(z)  -  0  (3.21) 

y  —  oo 

ILitil  we  reach  the  point  where  we  must  deal  with  the  radiation  condi¬ 
tion,  we  may  safely  omit  the  indices  c  and  s  for  the  calculations  are  the 
same  for  each. 

Let  us  now  introduce  the  auxiliary  function 

G(z)  -  f'(z)  4  ikf(z)  (3.22) 

which  according  to  (3.17)  satisfies  the  condition 

Im  G(x  4-  10)  -  0  (3.23) 
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By  the  Schwarz  reflection  principle,  we  extend  G  analytically  into  the 
upper  half -plane: 


G(x  +  iy)  -  G(x  -  iy) 


(3.24) 


or 

<t>x(x,y)  -  kH'Cx.y)  *  <|>x(x,-y)  -  k^x.-y) 
and  H>y(x,y)  +  k(j>(x,y)  «  <py(xt-y)  -  k$(x,-y) 

We  have  assumed  in  (3.20)  that  |f'|  <Bif  |*|  i  l  +  1,  y  <  0.  Let  us 
examine  what  this  may  imply  for  |g|.  Let  a  be  some  point  of  the  region 
(see  Figure  7). 


Figure  7  -  Path  of  Integration 


Then 


and  consequently 


f  (x) 


f ’ (*)dz  +  f (a) 


f ' (z)dz|  + 


!*(•>! 


<  B  |z  -  a|  +  j f (a) |  <  B  |z|  +  A 


But  than 


|G(x) |  <  C  |x|  +  D,  C,  D  >  0  (3.25) 

an  inequality  that  holds  in  the  whole  region  |x|  _>  A  +  1»  not  just  in  the 
lower  half -plane. 

Neer  the  point  x  •  -  iA,  we  have  assumed  (3.19).  Define 

M(x)  i  (*2  +  A2)1/2  f *  (*) 


Then 


M’(x)  "  2  *~2~T/2  f'  +  (*2  +  f" 

(s2  +  aV' 


or 


U2  +  A2)  M’  (*) 


s(s2  +  A2)1/2  f’  ♦  (*2  +  A2)3/2  f" 


Consequently 

<*2  +  A2)3/2  G’(z)  -  (*2  +  A2)3/2  f"  +  ik(s2  +  A2)3/2  f’ 

-  (s2  +  A2)  M'(*)  -  *  M(s)  +  ik(*2  +  A2)  M(s) 

=  N(x) 

Thus  (x2  +  A2)3^2  G* (x)  is  analytic  at  the  point  x  -  -iA  and  then  also 
through  the  reflection  at  x  ■  iA. 
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Sines  (s2  +  i2)3^2  G'(t)  is  snslytic  at  s  ■  +  if,  it  follows  thst 
2  2  1/2 

(s  +  t  )  G(s)  is  slso  snslytic  it  i  ■  +  iJt.  To  sss  this,  ons  may  start 
with  G* (s)  and  calculata 


G(*)  -  G(a)  -  J  (C2  +  i2)“3/2  H(C)  dC 


l2  (s2  +  i2)1/2 


N(«) 


i  x 

2  J  2  2  1/2  ^  ^ 

a  r  4  (C  +  * )  ' 


or 


(*2  +  i2)1/2  G<s)  -  (*2  +  i2)1/2  f*G(a)  -  ^  ~2  A  ~~2~i/2  N(a)l  +  'T  *N<*> 

L  i  (a^  +  in1'  J  r 

-  ij  («2  +  l2)1'2  J  (t2  +',2,1/2  (?)dt 


It  remains  to  be  shown  that  the  last  integral  is  bounded  near  x  ■  -11. 

This  is  uot  difficult,  and  we  shall  omit  the  proof  here. 

Let  us  now  develop  G(s)  in  a  Laurent  series.  Because  of  (3.25),  the 
series  will  take  the  fora 

G(s)  -  cs  +  aQ  +  ~  + . .  |  x |  >  l  +  1  (3.26) 


The  condition  (3.23)  implies  that  c  and  all  the  aQ  must  be  real.  Mow 


G*  (*) 


•  •  • 
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0  and  hence  that 


It  follows  aaslly  from  lia  f  »  0  that  lim  f"  ■ 

y  4  -  «o  y  4  -  n 

lim  G’(s)-0  (3.27) 

y  -4  —  oo 

But  than  in  (3.26) ,  tha  constant  c  -  0. 

Let  us  now  take  a  circle  C  with  center  a  ■  0  and  radius  greater  than 
1+1  (see  Figure  8) .  Hence  the  Laurent  series  is  valid  on  and  outside  the 


Figure  8  -  Region  for  Laurent  Series 


circle  and  the  Cauchy  theorem  yields 

4  G(t)dC  -  211^ 

■^c 

Since  G  is,  in  fact,  analytic  everywhere  outside  the  flat  plate,  the  inte¬ 
gral  along  C  can  be  deformed  onto  the  two  sides  of  the  plate  with  little 
circles  of  radius  e  around  the  ends,  where  G  may  have  square-root  singu¬ 
larities  according  to  what  we  have  shown  above.  Hence 

<»2-e  r  -ft+e 

2iria.  -  J  G(+0  +  in)  idn  +  J  G(-0  +  in)  idn 

-!-e+  £-e 
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+  I  G(-i£  +  Ee  )  Ee  id© 
-3n/2 


f3ir/2 

+  |  G(i£  +  ee10)  ei0id© 
-1/2* 


Since  G(+i£  +  eei0)  -  e-1^2,  the  contributions  from  the  small  circles 
vanish  as  £-+0.  We  have  from  the  definition  of  G 


f  +° 

2iria.  -  J  [A  -  Vi  -  i(4>  -  k*))  idn 

1  -£+€  x  y  -0 


14»  -  k<M  dn 


We  have  been  able  to  set  $x  -  0  from  (3,.  14).  Since  the  plate  must  be  a 
streamline,  it  follows  that  ¥(+0,n)  “  ^(-O.n)  *  C  and  hence  that  the  inte¬ 
gral  of  ¥  vanishes.  Since  the  left  side  is  purely  imaginary  and  the  right 
side  is  real,  we  can  conclude  at  once  that  both  must  be  zero.  However,  we 
can  also  show  that  the  right  side  is  zero  by  invoking  the  antisymmetry 
property  in  (3.24).  In  any  case  ve  have  shown  that  a^  -  0. 

Evidently 


2a,  3a. 

G’(z) - 3  4  “ 

z  z 


if  I  z  I  >  £  +  1 


(3.28) 


But  then  the  function  (z2  +  £2)^2  G*(z)  is  bounded  in  the  region  |*|  i  A  +  1 
We  have  also  seen  that  it  is  analytic  at  z  ■  +  i£  and  hence  also  in  the 
region  I  z  |  <  £  +  1  except  possibly  at  z  «  0,  where  there  is  a  confluence 


of  boundaries  with  different  boundary  conditions.  We  could  allow  a  singu¬ 
larity  at  z  ■  0  in  order  to  represent  a  loss  of  energy  in  breaking  of  waves 
on  the  barrier.  The  theory  would  not  tell  us  how  strong  to  make  the  singu¬ 
larity  but  it  would  serve  as  a  mathematical  model  for  an  energy  absorber. 

We  shall  not  do  this,  however,  but  instead  assume  that  there  is  no  loss  of 

2  2  3/2 

energy.  Consequently,  (z  +  i  )  G'(z)  is  bounded  in  the  whole  complex 

plane  and  hence,  by  the  Liouville  theorem,  must  be  a  constant: 

G'(z) - 5 - %  ■-/,  ,  c  real  (3.29) 

(z*  +  J/V' 

Integrating  once,  we  obtain 

G(z)  -  %  2  2  1/2  +  d*  d  rea3L  (3.30) 

£  (z  +  l  ) 

2  2  1/" 

Up  to  now  we  have  avoided  specifying  what  we  mean  by  (z  +  l  )  .  We 

choose  the  branch  that  behaves  like  z  at  large  distances  from  the  origin. 

2  2  1/2 

Then  on  the  right  side  of  the  flat  plate,  (z  +  i  )  takes  the  value 

(Z2  -  y2)1/2  and  on  the  left  side  -(l2  -  y2)1/2. 

There  will  be  no  loss  of  generality  if  the  constant  aQ  in  (3.26)  is 
chosen  to  be  zero  for  this  can  be  accomplished  by  adding  a  suitable  con¬ 
stant  to  f  itself  and  has  no  effect  upon  the  motion.  With  an  ■  0,  we  see 

-2  U 

from  (3.26)  that  G(z)  behaves  like  z  as  z  This  behavior  can  be 

9 

obtained  in  (3.30)  by  setting  d  ■  -c/Ji*’. 

Hence 

GU) ' c  ['(TT?)175 ' x] ' c  rMl  (3,31) 

This  completes  the  first  step  in  the  reduction  method,  the  finding  of  G. 

We  must  now  proceed  to  find  f. 
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Since 


G  -  f  *  +  ikf  -  e"ik*  [eik*  f<*>] 
we  nay  inucd lately  integrate  the  dif feren  ial  equation  for  f  and  obtain 


f(*) 


-ik* 

c  e 


TV75 


lj  «lkc 


In  order  to  be  specific  about  the  constant  of  integration,  we  shall  start 
the  integral  at  z  »  +  i<»  and  Integrate  along  a  path  that  lies  on  the  right 
of  the  flat  plate,  as  shown  in  Figure  9.  However,  in  order  not  to  lose 


Figure  9  -  Path  of  Integration  for  Equation  (3.32) 
generality,  we  must  now  add  a  constant,  i.e., 

f(*>  -  c  e'lk*  f  [  ~ *2 — T Tn  ~  «ikCdC  +  B  e~lk*  (3.32 

i»  lor  +  l  )  J 


In  order  to  determine  the  constant  B,  we  shall  use  condition  (3.15). 
For  this  purpose  we  oust  calculate 


f '  (z)  -  -ifcCe" 


J  [  2-^-2  172-4  *lk?d« 

i«  L(s  +  a  r'  J 


[  (z2  +  i2)1^2  x] 


Jz 

— dc  (3.3 

vz  ~r  *  ,  -ri*  (C“  +  tryu 


Then,  if  we  choose  a  value  of  y  on  the  right  side  of  the  plate,  we  obtain 


fU 

f’(+  0  ♦  iy)  -  -  ikBeky  +  C  ----  ** -  ikCe **  J  - T^-TTTJ 

(ir  -  y ) '  +i»  Kn  -  iV' 


-  ik.Ce 


ky  fiy  me“kn 
if  (i.2  -  n2)1 


-  ikn.ky  +  1C  y  y  TYfi  -  ko.ky  J  ;  t/2  dn 

{l*  -  yV'  l  (n  -  ir)l/ 


+  lkc.ky  C  — 


If  we  write  B  ■  B^  +  iBg,  then 

®  k 

Re  f * (40  +  ly)  -  kB_eky  -  kCeky  |  i  dn  -  0 

1  *1  (n  -  jrr' 

The  integral  is  a  known  one,  namely, 

[  (rf-l2)1'2  d"  “  1  Kl<k*') 

* 

where  is  a  modified  Bessel  function  in  Watson  notation.  But  then 
B2  “  C  %  K^(kf-)  and  we  have 

lj  eiktdC  +  tBx  +  iC*K1(kJ4)]  e~ikz 

(3.34) 

where  C  Is  real  (see  (3.31)).  This  function  now  satisfies  the  free-sujrface 
condition  (3.17),  the  condition  on  the  plate  (3.15),  the  sharp-corner  con¬ 
dition  (3.19),  and  the  conditions  at  infinity  (3.20)  and  (3.21).  We  must 
now  bring  into  che  picture  the  separate  functions  f and  fDg,  the  condi¬ 
tions  (3.16),  and  the  radiation  condition  (3.18).  However,  before  we  dp 
this,  it  will  be  useful  to  have  the  asymptotic  behavior  of  f(z)  as  x  -*■  +  00 . 


4 

See  G.K.  Watson,  "A  Treaties  on  the  Theory  of  Bessel  Functions," 
Second  Edition,  Cambridge,  at  the  University  Press  (1944) . 


For  this  purpose  we  deform  the  path  of  integration  as  shown  in 
Figure  10. 


Figure  10  -  Deformed  Path  of  Integration 


If  we  let  R  00  in  the  integral  along  the  quarter-circle,  it  is  not  dif- 
ficult  to  show  that  the  contribution  from  this  part  converges  to  zero. 
This  then  gives  us 


1  tcTT 


P~I72 


1 


elkCdc 


where  the  path  of  integration  in  the  second  integral  must  go  under  the 
plate.  If  we  now  let  x  -►  +  °°,  it  is  evident  that  this  integral  must  con¬ 
verge  to  zero  if  it  is  to  exist  as  at  improper  integral-  which  it  dors. 
Hence 


£(z)  ~  +  iC£K^(ki)J  e  as  x  •*  » 


(3.35) 


As  x  -  °o,  we  shall  have 
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f(z)  -  [B.  +  iCJIK  Od)J  e"ikz  -  C  e_ikz  J  -  1  eiki;dC 

11  [(5  +  *  ) 


where  the  path  of  Integration  is  below  the  plate  as  shown  in  Figure  11. 


Figure  11  -  Path  of  Integration  Below  the  Plate 


We  shall  complete  the  path  of  integration  by  adding  a  large  semicircle 
in  the  upper  half-plane,  as  shown  in  Figure  12.  Again,  it  is  easy  to  show 
that  as  R  +  ®,  the  contribution  from  the  integral  along  the  semicircle 


Figure  12  -  Completion  of  Path  of  Integration 


vanishes.  Hence,  adding  it  to  the  original  integral  has  not  changed  its 
value.  To  be  more  precise,  what  we  have  shown  is  that 
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where  CR  is  the  semicircle  of  radius  R.  Since  the  integrand  is  an  analytic 
function,  we  may  deform  the  integral  right  onto  the  extended  flat  plate  (see 
Figure  13).  Just  as  in  the  earlier  reasoning  concerning  G,  the  contributions 


1 

< 

> 

Figure  13  -  Path  of  Integration  Around  Extended  Flat  Plate 

from  the  little  circles  at  the  ends  go  to  zero  as  their  radii  go  to  zero, 
and  we  are  left  with 


Here  again  1^  is  a  modified  Bessel  function  of  well  known  characteristics. 
Finally,  then,  we  have 

f(z)  ~  [*1  +  iCSJ^ki)  -  2irACI1CkJl) ]  e  1KZ  as  x  +  -  «  (3.36) 
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We  arc  now  ready  to  invoke  the  radiation  condition.  Proa  now  on.  the 
function  f  and  the  conatanta  and  C  will  all  carry  additional  indicea 
c  or  a.  If  we  refer  to  (3.11)  and  (3.13),  we  aee  that 


Ic(.)  -  -i  £  .lkk  +  tx 


fe<*) 


M  e-ik*  +  f 

cr  DS 


Consequently,  if  we  apply  (3.35)  and  (3.36)  to  f  and  f  , 

C  8 


we  find 


fDC<*>  '  [1f4+»lc  +  *V,tl] 

£DS<*>  -  (f  +  «!.  +  £VKl]  «"ll“  ^ 

£dc‘*>  '  [£  f  +  Bu  +  lcctKi  -  “."l] 
fDS(l>  -  +  »1.  +  1Cc*Kl  '  2"C.lIl] 


88  X  +  +  » 


(3.37) 


Then  the  equation 


as  x  +  -  « 


lim  (f ’  +  kf__)  -  0 

x  S 


(see  (3.13))  gives 


-  lk  [£  f4  +  *lc  +  1CcUl]  +  k  [S4  +  »1.  +  £C.“l] 


and  the  equation 


gives 


,  il!  „(tDC  •  k,DS>  ■ 0 


-  lk  [‘  £  *  llc  +  lCc“l  -  *“c«l]  -  k  [f  +  *U  +  lC.lKl  -  2,c.lll]  •  0 

Separating  real  and  imaginary  parts,  we  find  the  following  set  of  four 
equations  for  the  unknown  constants  B^fi,  ,  C^c>  C^: 


-B 


lc 


•  v 

1  c 


-  ULC#  -  0 


-B 


la  +  JU^C,  +  2x11^  -  0 


B 


lc 


-«ic. 


B 


is 


+  1K.C 
1  c 


_  2gA 

o 


The  solution  of  these  equations  yields 

St 


Kl(kA) 

Cc  *-  7  2  .  -2,  2  at  * 


lc 


K1  +  1Tll 


~7rIl>Cl  jgA 

Kl2+irV  0  ’ 


c  -  - 


*1  !<W)  ^ 

Kl2  +  7f2ll2 


'Is 


■  l/:  »v  • 2]  * 


Substituting  into  the  asymptotic  expressions  (3.37),  we  find 


DC 


.  **  J  i-i _ +  i -lh _ 1 

LKi  +ui  Ki  +  niJ 


1#  m 

i  a»x  +  +  «  (3.38) 
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f  -  ^ 
dS  o 


r  ,y  w.rt  I 

U2-V  V-vJ 


e“ik*,  as  x  -*•  +  • 


f  «. 

DC  O 


f  _ 

DS  o 


r  ni*i  , .  *v  i  -it, 

L^V  ^V J' 

r  i  ri  . ,  i.-u« 

W + *v  \2 + *  v  J 


The  sum  f  "  fDC  +  fIC  and  fTS  ”  fDS  +  flS  *ive  the  two  parts  of  the 
transmitted  wave  when  x  <  0.  The  asymptotic  expressions  are  as  f ollct’s: 


m  [—hh~ .  t  ..  ’X.  .  1 
LA  +  *1  *1  +  ”  J 

.  it  r__aL. .  t ,  whh .  I.- 
°  U2+’V 


as  x  +  • 


(3.39) 


Let  us  use  the  results  shown  to  compute  tj  u  x  +  +  »,  the  velocity 
potential  of  the  reflected  wave,  and  $T  «  ♦j  +  ♦p  as  x  -  ®,  that  of  the 
transmitted  wave: 

♦D  -  eky  — j — ^~2 — 2  1TllKl  co*  +  ***1*  ain  cos  at 


+  vl 


2  2 

+  [—  ti  1^  coo  kx  -  sin  kx]  sin  at 


(3.40) 


.M.w 

a  I  2  .  2.  2 

lh  +vh 


it\2  J 

cos  (kx  -  ot)  +  — = - = — j  sin  (kx  -  at)  I 

K1  +  "  V  -J 

II  x  +  +  * 
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1A  ky  I"  **1*1 _ 

W+*i 


2  co*  (kx  +  at)  -  — j - 2 — 2  *in  +  ot) 


K:  +  ir  lx 


as  x  - 


(3.41) 


Th«  associated  expressions  fer  the  free  surface  are 


yd- A 


rA!_ 

LKi2  +  *V 

[7^ 


vhh 

cos  (kx  -  at)  +  — 2 - 2 — 2  *in  ^  “  ot^  1  »  x  +  ** 

K1  ♦  *  xi  J 

(3.42) 


K.  Tt  K. 

YT  ~  A  — 2 - 2 — 2  cos  +  0t^  +  — 2 — ~2 — 2  8in  +  at>  ,  x  » 

+  n  Kj  +  tr 


It  Is  customary  to  Introduce  a  reflection  coefficient  R  -  (Amplitude 
of  Yd) / (Amplitude  of  Yj)  and  a  transmission  coefficient  T  -  (Amplitude  of 
Y,p) / (Amplitude  of  Yj),  Since  here  (see  (3.9))  the  amplitude  of  Yj  -  A, 
we  find 


^  +  *V>1/2 


^(tt) 


T  "  2  .  072  (3*43) 


QLf  +  *  Ip 


The  equation  F  +  -  1  is  simply  an  expression  of  the  conservation  of 

energy.  If  we  had  allowed  a  singularity  at  s  -  0  to  represent  wave  break¬ 
ing,  this  equation  would  no  longer  hold. 

We  define  the  phase  shifts  by  comparing  the  actual  reflected  wave  with 
the  completely  reflected  wave  A  cos  (kx  -  at)  and  the  transmitted  wave  with 
the  incident  wave  A  cos  (kx  +  at): 


*  -*■  ♦ 


*D  ~  At  cos  Out  -  ot  -  , 

YT  ~  At  COS  (k*  01  -  eT) 

Ms  easily  find 

K-  TTl. 

c“  •  tM  st  *  (J- 

Consequently 

S  +  PT  "  7  n  (3. 


The  functions  and  1^  are  wall  tabulated,  so  that  there  is  no  dlf 
flculty  in  plotting  A,  T,  and  ^  as  functions  of  kl.  The  graphs  look 
approxlnstely  as  shown  in  Figure  1A.  It  Is  evident  that  if  X  «  l,  the 
reflection  is  * 


* 

* 


almost  complete  vinrui  if  \  »  l,  the  transmission  is  alaost  coaplttc, 
rasults  that  accord  with  our  intuition.  However,  comparison  with  experi¬ 
mental  measurements  doas  not  show  vary  good  agreement.  Part  of  tha  raason 
for  this  is  probably  naglact  of  tha  vorticas  that  ara  fornad  by  tha  flow 
back  and  forth  around  tha  sharp  adga  at  tha  bottom  of  tha  plats.  Thair 
formation  is  not  takan  into  account  in  tha  formulation  of  tha  problem. 

Wa  hava  baan  discussing  only  tha  asymptotic  behavior  of  tha  solution 

whan  x  -*■  +  •.  Hovevar ,  as  soon  as  wa  hava  found  Cc,  and  C^,  wa 

nay  substitute  into  (3.34)  and  hava  tha  complete  solution  for  f#(x)  and 

f  (s).  With  these  functions,  wa  can  than  construct  tha  solution 
c  io 

for  tha  whole  region  y  <_  0.  In  particular,  wa  can  find  the  pressure  on  the 
two  sides  of  the  plate  and  calculate  the  force  and  moment  about  0  acting  on 
the  plate.  We  shall  not  do  this  to  avoid  getting  bogged  dawn  in  treating 
just  one  problem.  However,  this  has  been  done  by  Haskind  (1948);  a  brief 
discussion  of  the  results  is  given  in  We hausen  and  Laitone  (1960,  pp.  532- 
533). 
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PART  4:  THE  METHOD  OF  GREEN  FUNCTIONS 

Of  the  several  methods  that  ve  ere  discussing  here,  the  method  of 
Green  functions  is  certainly  the  most  flexible  in  application,  allowing 
treatment  of  a  much  wider  class  of  problems  than  the  other  methods. 

However,  even  though  it  has  been  known  for  a  long  time  and  has  given  rise 
to  many  mathematical  investigations,  it  did  not  really  become  important  in 
the  solution  of  engineering  problems  until  perhaps  the  last  ten  years.  The 
reason  for  this  is  that  a  "solution"  by  this  method  typically  involved  the 
solution  of  an  integral  equation.  Although  one  could  show  in  msny  cases 
that  a  solution  to  the  integral  equation  existed,  an  analytic  solution  was 
usually  not  obtainable  and  a  numerical  solution  was  too  difficult.  The  ad¬ 
vent  of  high-speed  computers  has  radically  changed  the  situation  with 
respect  to  numerical  solutions,  and  nowadays  this  is  almost  routine. 

The  above  remarks  should  not  be  Interpreted  as  meaning  that  one  cannot 
obtain  an  explicit  solution  by  using  Green  functions.  However,  in  those 
situation?,  where  one  can  do  this  one  can,  as  far  as  1  am  aware,  also  solve 
the  problems  by  another  method.  For  example,  Havelock  derived  the  Michell 
integral  by  using  a  Green  function.  However,  Michell  himself  derived  it 
from  Fourier  analysis,  having  first  derived  a  representation  theorem 
analogous  to  (2.22). 

We  8 ha 11  approach  this  method  by  way  of  some  Green  identities,  which, 
in  turn,  are  directly  derivable  from  the  Gauss  divergence  theorem.  Let  <$> 
and  4*  be  any  two  functions  defined  in  a  certain  three-dimensional  region  V 
and  having  second  derivatives  there,  including  t;®' boundary  S.  Then  one  of 
the  Green  identities  states  the  following: 

dV  -  I  [W  -  W  1  dS  (4.1) 

Js  n  n 

where  the  normal  vector  points  out  of  the  region  V.  If  in  addition  both  4> 
and  <|>  satisfy  the  Laplace  equation,  the  left-hand  side  vanishes  and 


I 


IW_  -  <M»  ]  dS 
n  n 


(4.2) 
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A  particular  solution  of  the  Laplace  equation  in  three  dimensions,  in 
either  P  -  (x,y,z)  or  Q  -  (C.h.O,  is 

1/2 

¥  -  7  ,  r  -  [(x-o2  +  (y-n)2  +  (z-O2]  (4.3) 

The  only  exception  is  at  the  point  P  ■  Q,  where  the  solution  is  singular. 

Let  us  take  Q  as  the  integration  variable.  We  must  then  exclude  a  small 
region  containing  P  from  the  integration.  Let  this  be  a  sphere  V£  of  radius 
£  with  surface  S£.  Then  (4.1)  takes  the  form 


Since  the  left-hand  side  vanishes  identically  and  since  the  integral  over 
S£  converges  to  -4  it  <J>(P)  as  e  0,  we  find  another  identity: 


(4.4) 


Here  we  have  written  v  instead  of  n  in  order  to  make  clear  that  the  normal 
derivative  is  taken  in  the  variables  (£,m,4). 

If  the  point  P  is  taken  on  the  surface  S  at  a  point  where  the  surface 
is  smooth,  we  may  show  that  (4.4)  stili  holds  if  4tt  is  replaced  by  2it. 

We  may  now  generalize  (4.4)  further  by  exploiting  (4.2).  If  we  add  to 
1/r  any  harmonic  function  y  (i.e.,  one  satisfying  the  Laplace  equation  and 
having  no  singularities  in  V),  (4.4)  will  still  hold,  We  shall  choose  a 
special  kind  of  function  to  add  to  1/r.  Let  H(x,y,z,C,T|,^)  be  harmonic  in 
each  set  of  variables  (x,y„z)  and  (£,n,C)»  We  define 
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G(P,Q)  -  ±  +  H(P,Q) 


’n.-gj1 1  wwp 


(4.5) 


1 


* 


It  then  follows  from  (4.4)  and  (4.2)  that 


J 

'■I 


4>(P)  “  ^  f  [^<Q)G(P,Q)  -  <j>G^]dS(Q) 
Js 


(4.6) 


As  before,  if  the  point  P  is  taken  on  S,  the  4ir  is  replaced  by  2ff. 

We  have  been  considering  only  three-dimensional  motion.  However,  there 
are  corresponding  theorems  if  the  motion  is  two  dimensional.  In  this  case 
(4.5)  is  replaced  by 


G(P,Q)  =  log  i  +  H(P,Q) 


(4.7) 


■ 


2  2 

where  P  =  (x,y),  Q  *  (£,ri)  and  r  ■  [(x-£)  +  (y-H)  )  •  Formula  (4.6) 

becomes 


r 


4>(p)  m  h  f  [V  ■  *Gv]ds 

Jc 


(4.8) 


where  C  is  a  contour  bounding  a  two-dimensional  region  S  and  dd  is  arc 
length.  If  P  is  on  C  at  a  point  where  C  is  smooth,  then  the  2tt  is  replaced 
by  tt. 

Let  us  now  turn  to  several  problems  where  we  can  exploit  (4.6)  or  (4.8). 


FORCED  HARMONIC  MOTION 

Let  us  suppose  that  a  body  is  being  forced  to  undergo  periodic  motion 
of  frequency  a  and  in  such  a  way  that  there  is  no  average  displacement. 

We  shall  denote  the  average  position  by  Sq.  Let  the  Oy-axis  pass  through 
(or  near)  the  body  and  let  Z  be  a  cylinder  of  radius  R,  large  enough  to 
contain  the  body.  We  further  denote  the  portion  of  the  free  surface  F 


fc  -Li  :7V*- 
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inside  XR  by  Fr  and  the  portion  of  the  bottom  ft  inside  i.R  by  (sea 
Figure  IS). 


If  we  write 

<&(x,y,z,t)  »  ^(x.y.z)  cos  at  +  sin  at  »  Re  $e  iat  (4.9) 

where  $  »  then  from  (1.13),  the  free-surfece  boundary  condition 

4>y(x,0,z)  -  \xj>  -  0,  v  ■  a2/g  (4.10) 

The  condition  on  the  bottom  is 


and  the  condition  on  Sq  is 


0 


-  V 

n 


(4.11) 


(4.12) 


Since  this  is  still  in  terms  of  $  rather  than  4>,  let  us  decompose  V  : 


V  (x,y,z,t)  -  V  .  cos  Ot  +  V  ~  sin  ct 


Ra  Vo(x,y,*)  s 


(A.  13) 


where  V  (x.y.s)  »  V  ,  >  IV  Usually  the  face  that  the  body  dlaplacaaents 
n  nl  nz 

are  saall  1c  exploited  further  and  we  can  write 


Vo(xty,x,t)  -  ^(t)  nx  ♦  «*3 


♦  Afrxn)^  +  8(rxn>2  f(rxn)j 


(A. 14) 


Here  (x^(t) ,y^,z^)  describe  the  translation  displacements, 

(u(t),8,Y)  describe  the  angular  diaplac— enta  of  the  body,  and  r  is  a 
vector  fro*  the  origin  to  a  print  (x,y,z)  on  the  body  surface  S„.  If 
we  define 

ai  "  V  a2  •  yr  a3  •  *i 


-  0z  -  yy,  a5  -  yx  -  a*.  *  <*y  -  Sx 


we  nay  write  (4.14)  in  the  fora 


(4.15) 


Vn(x,y,z,t)  -  y  (^(t)  ^(x.y.zj 


(4.16) 


Now  afe(t)  itaeif  can  be  written  as 


ofe(t)  1  a^'  cos  ot  +  a^"  sin  ot  »  Re  *> 


(4.17) 


The  boundary  condition  (4.12)  now  gives 


(k)  |  '  *  00 , 

n  U  nV»  » ! 


(4.19) 


If  we  wish  to  divide  <t  into  sine  and  cosine  components,  we  may  write  it  as 
follows : 


*  "  ZW  0lCk>  +  V  4,2(k>3cOS  qt  +  y.al-  V  ^l(k)  +  V  <J>2(k>lsln  0t 

k  k 

so  that  in  the  notation  of  (4.9) 

*1  -  Z7atak"  ♦lW  +  V  *2 01)3 

k  (4.20) 

*2  "  *l<k)+\^2(k)> 


To  justify  the  decomposition  (4.18),  we  should  demonstrate  that  the 
problem  has  a  unique  solution.  It  is  possible  to  do  this  after  we  add  one 
more  condition,  the  radiation  condition.  In  the  two  examples  treated  by 
the  Fourier  method,  we  satisfied  this  by  manipulating  some  undetermined 
constaata.  However,  when  treating  the  diffraction  about  a  plate  by  the 
reduction  method,  we  formulated  the  radiation  condition  directly  in  terms 
of  f^,  and  f0J.  (see  (3.18)).  There  ia  an  analogous  formulation  for  the 

present  case.  If,  analogously  to  (4.9),  we  write  <^k3  -  <$>.  ^k3  +  id./k\ 

(k)  1  ^ 

then  the  radiation  condition  for  '  takes  the  following  form: 

lim  *l/k  Id  00  -  ik  d(k)l  -  0, 

0  2  (4.21) 

<r  »  gkg  tanh  k^h 

2  2 

where  K  *  (x  ♦  a  J  .  The  radiation  condition  for  ♦  «  ^  +  i$2  !•  the 
aame.  It  ia  a  consequence  of  (4.21),  which  wa  shall  not  prove  here,  that 
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(4.22) 


'  -  0(R  '  )  as  R  -►  » 

i.e.,  that  <M<“»ifR>RQ  for  some  sufficiently  large  Rq. 

As  far  as  boundary  conditions  are  concerned,  the  only  difference  be- 
(k) 

tween  4>  and  the  4>  is  in  the  boundary  condition  on  Sq.  All  of  thee 
satisfy  (4.10),  (A. 11),  and  (4.21).  However,  satisfies  (4.19)  whereas 

<j>  satisfies 

I e  ■  V  (x,y,*)  (4.23) 

n  j>0  n 

Much  of  the  following  can  apply  to  either  one.  We  shall  write  out  the 

(k) 

development  for  <j>.  To  convert  to  4>  ,  we  need  only  replace  Vn(x,y,e)  by 

nk(x,y,z). 

Let  us  now  apply  (4.6)  to  the  domain  bounded  by  Sq,  Fr,  Br,  and  ZR. 

We  find  the  following: 


4ir<j)(P) 


1 


ouFruBruEr 


14>V(Q)  G(P,Q)  -  <KQ)  Gv(p,Q)]dS(Q) 


[Vv(Q)  G(P,Q)  -  <KQ)  GJ(P,Q)]dS 

2 

[j~  <KC,0,O  G(P;C,0,C)  -  <J>  Gn  ] 


-  $(Q)  Gv(P,Q)  dS 


[$rg  -  <j)GR]  Rdedn 


(4.24) 


where  we  have  used  (4.23)  in  the  first  integral,  (4.10)  in  the  second,  and 
(4.11)  in  the  third. 


We  now  consider  several  cases 


Casa  1.  Ve  begin  by  writing  the  last  integral  of  (4.24)  as  follows: 


1 


(R1/2  [*R  -  ikQ^l  R1/2G  -  R1/2  [Gr  -  ikQG]  R1/24»}  d8dn 


(4.25) 


“R 


,1/2, 


In  this  Integral,  we  know  from  (4.21)  and  (4.22)  that  R  [<j>_  -  1K$)  0 

1/2  K 
and  that  R  $  rewains  bounded  as  R  «.  Let  us  now  suppose  that  it  is 

possible  in  (4.5)  to  find  a  function  H  such  that 


Gn(P;S.o,c)  -  \»G  -  0 


Gv(P,Q)  •  0  for  Q  on  B 


(4.26) 


1/2  2  2 
lie  R1  ^[GR(P»C»n»C)  -  ik0G]  -  0,  R  -  a  +  C) 


1/2 


Any  consequences  we  may  draw  frow  these  properties  depend,  of  course,  on 

our  being  able  to  construct  such  a  function. 

It  is  now  immediately  evident  that  the  integrals  over  FR  and  BR  ere 

sero.  In  the  last  integral,  in  the  form  (4.25),  we  see  that  now 
1/2  1/2 

R  (GR  -  ik,jG]  0  and  R  G  remains  bounded  as  R  -*■  •.  Consequently  the 
Integral  over  ER  converges  to  sero  as  R  -*■  ®.  This  leaves  us  then  with  only 
the  integral  over  Sq.  If  we  now  let  P  be  a  point  of  Sq  instead  of  a  point 
in  the  interior,  the  4n  is  replaced  by  2tt  and  ve  nay  write  the  resulting 
equation  in  the  following  fora: 


♦  <P)  +  H  J  <KQ)Gv(P,Q)dS(Q)  -  ^  J  Vv(Q)G(P,Q)dS(Q) 

so  so 


(4.27) 


This  is  an  integral  equation  for  the  determination  of  $(P)  for  PeSg.  If 
this  equation  has  been  solved,  then  (4.24),  which  takes  the  form  below, 
gives  $(P)  for  any  point  inside  the  fluid: 
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where  and  X^  are  often  respectively  called  "added  masses"  and  "damping 

coefficients."  In  this  case  they  are  functions  of  the  frequency  as  well  as 

the  body  geometry.  There  are  some  advantages  for  further  developments  in 

(k) 

replacing  n^  in  (4.30)  by  ,  according  to  (4.50).  However,  we  shall 

net  explore  this  further. 

Under  what  conditions  can  one  find  a  Green  function  that  satisfies 
(4.26)?  It  is  possible  to  construct  such  a  function  for  either  infinite  or 
constant  finite  depth  by  applying  either  the  Fourier  method  or  the 
reduction  method.  The  functions  can  be  found  in  several  places,  e.g.,  in 
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If  P  Is  taken  near  Q,  it  is  evident  that  the  denominator  becomes  small  and 
1/2 

hence  that  R  It  becomes  large.  Properly  we  should  examine  the  integral 


whsru  K(p)  is  the  complete  elliptic  integrel.  We  can  undoubtedly  do 
better  then  this.  In  any  case,  it  is  not  enough,  for  K(p)  ■*  00  as  p  ♦  1. 

We  shell  approach  the  treatment  of  the  integrel  over  IR  somewhat  differently 
In  discussing  the  Fourier  method,  we  found  solutions  of  the  form 
(2.17).  If  in  the  solution  for  xQ  we  had  Introduced  polar  coordinatea 

x  ■  R  cos  0  ,  z  ■  R  sin  6  (4.35) 

we  should  have  found  the  following  solutions: 


Jn<k0R)] 

1 

>  x  < 

|  cos  n6 

j 

(sin  n6 

J  and  Y 
n  n 

are  Bessel  fi 

i 

(4.36) 


tation.  If  we  consider  the  region  exterior  to  some  bounded  body,  we  find 
that  the  only  contributions  for  x^  are  exponentially  decreasing  in  R  and 
that  the  part  of  a  solution  representing  outgoing  waves  has  the  following 
form: 


<j>  ~  a„  cosh[kQ(y  +  b)]  cos(n0  +  6r)  lJn(k0R)  +  iYn^ 

n»0 


where,  as  usual,  +  iY^ 


,(D 


.  The  combination  is  usually  called 


the  Hankel  function  of  the  first  kind.  Its  asymptotic  expansion  for  large 
k  R  is  well  known,  and  we  may  write  the  asymptotic  expression  for  $  as 


0 

follows: 


X 


♦  - 


^  «n  co«h[kg(y  ♦  h)l  coo(n6  +  6q) 
n-0 


(4.3R) 


|«»  -  f 


k  ntraight.forward  computation  yield* 


♦r  '  H  +  lk0*  ♦  '>Uk0R)*5/2l  (4.39) 

We  now  substitute  this  expression  for  into  the  integral  over  1  in 
its  original  fora  in  (4.24),  neglecting  ths  tera  of  0{OiqR)  )•  After 
soae  reordering  of  (4.24),  we  can  write  it  in  the  following  fora: 

dS 

(4.40) 

I*  i]  « 

•  I  Vv(Q)G(P,Q)dS 

''s 

b0 

If  we  now  let  P  approach  any  boundary,  we  have  the  ease  equation  with  4ir 
replaced  by  2n.  This  1*  then  an  integral  equation  for  the  determination 
of  4{P)  for  P  on  j£^.  Once  it  has  been  determined,  the  equation 

above  determines  $  at  any  interior  point.  We  should  keep  in  mind  that  in 
deriving  Fquation  (4.40),  an  approximation  has  been  made  In  evaluating  the 
integral  ova.'  1^. 

This  method  has  been  used  for  numerical  calculations  in  both  two  and 
three  dimensions  by  R.W.-C.  Yeung  (1973;  see  also  Bat  and  Yeung,  1974).  In 


4twKP)  +J  MQ)Gv(P,Q)dS  +J  MQ)  [cn(P»Q.'  -  -  o] 


►J  ♦(Q)Gv(P.^)dS  CQ>  [(iko’lg)  ~  ' 
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comp  wing  this  procedure  with  the  earlier  one  la  which  we  used  e  Green 
function  satisfying  (4.2ft),  the  following  points  should  be  kept  In  Bind: 

1.  The  Green  function  satisfying  (4.26)  la  such  store  complicated  than 
(4.31).  However c  this  Is  offset  by  the  fact  that  4  is  solved  only  for  P  on 
Sq.  In  the  Yeung  nethod  one  finds  6  for  P  on  SQuFRuBRuIR,  which  means,  of 
course,  a  much  larger  matrix  in  the  discrete  version  of  the  integral 
equation. 

2.  Although  in  the  Yeung  nethod  we  have  required  a  flat  bottom  for  R  large 
enough,  it  does  not  need  to  be  flat  in  the  region  near  the  body.  Further- 
more,  in  the  two-dimensional  version,  the  flat  region  can  have  different 
depths  on  the  two  sides.  This  did  not  seem  to  be  possible  with  the  more 
complicated  Green  function  satisfying  (4.2G).  We  now  turn  our  attention 
to  this  question. 

Case  III.  He  now  suppose  that  our  Green  function  satisfies  (4.26)  except 
that  G^-0  for  Q  on  B  is  replaced  by 

Gn(P;C,  -  h,e>  -  0  (4.41) 

We  have  mentioned  earlier  that  it  is  possible  to  construct  such  a  Green 
function.  This  Green  function  will  now  be  applied  to  the  situation  indi¬ 
cated  schematically  in  Figure  15.  However,  we  shall  suppose  that  the  part 
of  the  bottom  that  la  not  at  depth  h  is  limited  to  a  finite  stretch  that  we 
shall  denote  by  Bq  (see  Figure  16). 


J  > 
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If  we  now  use  this  Green  function  in  (4.24)  end  take  the  limit  R  -►  «, 
we  obtain  the  following  equation: 


4if$(P) 


-J  (Vv(Q)G(P,Q)  -  *(Q)GV]  dS 
JS 


+  f  -  4><Q)GV(P,Q)  dS 
B0 


(4.42) 


If  we  now  let  P  be  a  point  of  either  Sq  or  Bq,  we  obtain  the  following  in¬ 
tegral  equation: 


<!>(P) 


4>(Q)GV(P,Q) 


dS 


S0uB0 


Vv(Q)G(P,Q)  dS 


(4.43) 


This  is  almost  exactly  the  same  as  the  integral  equation  (4.26)  except  that 
the  integral  on  the  left-hand  side  is  now  over  SquBq  instead  of  just  SQ. 

Now  by  comparing  Cases  1,  II,  and  III  when  the  bottom  is  uneven  near 
the  body,  we  may  make  the  following  observation.  Case  I  is  not  feasible 
for  we  do  not  know  how  to  construct  the  Green  function.  Both  Cases  II  and 
III  are  feasible.  Case  II  has  a  simple  Green  function,  but  we  must  find 
<j>(P)  over  an  extended  boundary.  Case  III  has  a  more  complicated  Green 
function,  but  we  must  solve,  for  <£(?)  over  a  more  restricted  boundary. 

Some  comparisons  of  computer  times  for  the  case  of  a  completely  flat 
bottom  are  given  in  a  paper  by  Bai  and  Yeung  (1974) .  In  the  case  con¬ 
sidered  there,  a  two-dimensional  one,  about  twice  as  much  time  was  required 
for  tne  Yeung  method  (Case  II)  as  for  the  Frank  method  (Case  I). 


DIFFRACTION  PROBLEMS  AND  THE  HASKIND 
RELATIONS 

The  following  remarks  are  a  sort  of  appendix  to  Case  III  discussed 
above.  With  the  same  geometry  considered  there,  let  us  suppose  that  the 
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body  is  fixed  but  that  there  it  e  known  incident  wave  with  velocity  po¬ 
tential  ♦  .  If  ♦  ia  the  velocity  potential  for  the  total  fluid  notion, 
then  the  diffracted  weve  ip  defined  «e  usual  by 


(4.44) 


The  diffraction  potential  must,  of  course,  satisfy  the  free-aurfece  con¬ 
dition  and  radiation  conditions.  In  addition,  it  eust  satisfy  the  boundary 
condition 


Dn 


t>  on  both  Sn  and  B_ 
In  0  0 


(4.45) 


The  force  acting  upon  the  body  as  a  result  of  the  presence  of  the  fluid 
motion  is  given  by 

Fi  *  j  P“ids  “  “  P  J  (♦u  ♦  ♦Dt>nidS 

S°  S°  (4.46) 

-Re  |  -  iap  J  (^  -I-  $D)njdS  e”l0t j 


The  Haskind  relations  as  usually  presented  for  a  horisontal  bottom  or  in¬ 
finite  depth,  allow  solution  of  the  diffraction  problem  to  be  avoided  if 
the  forced-notion  problen  has  already  been  solved.  Here  we  shall  show  that 
this  is  also  true  if  the  bottom  is  uneven. 

Instead  of  the  function  $  defined  by  the  integral  equation  (4.43),  we 

(k) 

shall  use  the  function  '  defined  by  the  same  equation  with  V  (Q) 

*  (k)  v 

replaced  by  n^(Q).  Then  satisfies  the  boundary  condition  (4.19)  and 

in  addition 


(>.  -  0  on  B 

In 


(4.47) 


We  suppose  that  this  function  has  been  found. 
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i  J^r?***! i^rsjwqr -? m,  i*  !,i  i 1 


Consider  now  the  volume  bounded  by  S.UF  UB  UE  .  We  apply  the 

...  u  R  R  R 

Green  identity  (4.2)  to  4>pt  and  4>^  ' ; 


‘Wlf’  -  V^'l  dS 


SOuFRUBRu£R 


(4.48) 


[A  A  ^k)  a  <t>(k)  i  ds 

l^D9ln  9Dn?l  J  a 


soubouEr 


(k) 

for  both  <J>  and  <})1  v  equal  zero  on  the  flat  part  of  B  (recall  that 
*  0  there).  Also,  <j>^  -  0  on  B^.  Since  both  <J>^  and  <p^k;  satisfy 

the  radiation  condition,  the  integral  over  £R  tends  to  zero  as  R  ■*■  ®. 
Hence  we  may  write  the  last  equation  in  the  following  form: 


J  Ws  ■  J  Vi^k)dS  ■  J  v{k>ds 


S0UB0 


(4.49) 


‘•J 


4,In*(k)ds 


souBo 


Then  Equation  (4.46)  may  be  rewritten  as: 


F±  -  Re{-  lop  [  j  (*1»1<k)  -  *In*‘k))1S  -  J  4>In*<k,ds]  e10t) 


(4.50) 


is  now  expressed  completely  in  terms  of  known  functions.  The  ex¬ 
pression  (4.50)  differs  from  the  usual  one  only  in  the  presence  of  the 
integral  over  B^.  This  poses  no  problem  in  principle  for  in  solving  (4.43), 


we  have  found  '  on  both  Sq  and  B. 


3 


DIFFRACTION  ABOUT  ISLANDS,  HARBOR 
OSCILLATION 


As  a  sort  of  diversion,  I  should  like  to  consider  a  slightly  different 
type  of  Green-function  problem,  one  that  is  two  dimensional  in  the  mathe¬ 
matical  sense.  We  begin  by  considering  the  physical  situation  shown  in 
Figure  17: 


Figure  17  -  Representation  of  the  Physical  Problem  of 
Diffraction  about  Islands 


The  bottom  is  assumed  to  be  horizontal  and  the  island  to  have  vertical 
walls.  As  a  consequence  of  this  simplified  geometry,  the  problem  may  be 
treated  as  two  dimensional.  The  incoming  wave  will  have  the  velocity 
potential 


^jCP.t)  -  q  cosh'  k  h cosh  ko^z  +  sin  ^k(/x  cos  &  +  y sin  ^  ~  ot] 

°  (4.51) 

-  Re  4>I(P)e“iCTt  ,  O2  -  gkQ  tanh  kQh 


and  the  surface  profile 


Yj(x,y,t)  *  A  cos[kg(x  cos  3  +  y  sin  8)  -  ot]  (4.52) 

Note  that  here  we  have  taken  z  vertical  (see  Figure  18) . 

z 


x 


Figure  18  -  References  Axes  for  Diffraction  about  Islands 


The  diffraction  potential  ^(P.t)  *  Re  <J>g(P)  e  *at  will  satisfy  the  free- 
surface  condition,  the  radiation  condition,  the  bottom  condition,  and 


This  is  known  as  the  Helmholtz  equation.  Henceforth  we  shall  drop  the 
circumflexes 

As  a  preliminary  step,  we  go  back  to  one  of  the  Green  identities  in 
two  dimensions. 

J(uAv  -  vAu)dS  -  J  (uvn  -  vun)ds 
S  C 

If  both  u  and  v  satisfy  the  Laplace  equation,  the  left-hand  side  obviously 
vanishes.  However,  if  both  u  and  v  satisfy  the  Helmholtz  equation,  the 
left-hand  side  also  vanishes.  Hence 


'CT'?r^'TrTrw^*;<T,jvT»  v+'fr 


[Tyj.’.  '.1 1"'  V1  •“TfJHP r 


We  return  now  to  the  Island  problem.  Consider  the  region  bounded  by 
C  and  a  large  circle  of  radius  R,ER  (see  Figure  20).  We  hope  that  we  can 
construct  a  function 


/ 


sp 


\ 


\ 


Q 


\ 

Y 

L, 


\ 


/ 


Figure  20  -  Region  Bounded  by  C  and  £R 


G(P,Q)  -  log  ±  +  F 


-  log  i  +  Fx  +  1F2 

such  that  G  satisfies  both  the  Helmholtz  equation  and  the  radiation 
condition 


(4.57) 


/c  (G  -  iK.G)  -  0 
r-*»  r  0 


(4.58) 


We  shall  write,  according  to  the  formula  we  have  just  derived. 


Vp>  •  if  f  '*DVC  -  *DGVldS  +  if  /  l*0,C  -  Wds 

r  y. 


Since  both  and  G  satisfy  the  radiation  condition,  we  can  show  as  we 


did  earlier  that  lim  |  ■  0.  Hence 

R-*» 


I 


JR 


VP>  ‘ 


j(V>(Q)G(P’Q)  -  VvldS 
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or,  using  the  boundary  condition  on  C, 


VP)  +Tn  f  VQ)Gv<P,Q)dS  “  J  ^IV(Q)G(p.Q)dS  (A. 59) 

JC  Jc 

If  the  point  P  is  taken  on  the  boundary  C,  then,  as  before,  the  2tt  is 
replaced  by  it: 

VP)  +n  f  VQ)Gv(P,Q)dS  "  "  ¥  J  <J>IV(Q)G(P,Q)dS  (4.60) 

C  C 

This  is  an  integral  equation  for  ^  and  the  remaining  problem,  aside  from 
numerical  computation,  is  to  determine  whether  the  required  function  G  can 
be  constructed. 

We  shall  simply  give  the  function  rather  than  actually  construct  it. 

It  is  possible  to  show  that 

G(P,Q)  -  i  ~  [J0(kQr)  +  iY0(kQr)]  -  i  f  H<1} (kQr)  (4.61) 

The  integral  equation  for  is  then 

Vp>  +  i  J  VQ)Hov1)ds  -  -  i  J  *ivH™  (kor)dS  (4-62) 

c  c 

Once  <j>D  has  been  found  on  the  boundary  C,  then  we  can  use  the  more  general 
equation  (4.59)  to  find  ^(P)  for  P  outside  C: 

VP)  “  ~  i  J  ^uCQ)^^1^^)  +  <|.Iv(Q)H^1)(k0r)]dS  (4.63) 

C 

The  above  results  are  all  well  known  in  the  theory  of  diffraction  of 
acoustic  and  electromagnetic  waves  upon  infinitely  long  cylinders.. 

Let  us  now  consider  a  slightly  different  geometrical  configuration. 
Figure  21  is  supposed  to  represent  a  coastline  with  a  harbor  and  an  incident 
wave.  We  again  suppose  the  bottom  to  be  at  depth  h  and  all  coastlines  to 
be  vertical  cliffs. 
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Figure  21  -  Incident  Wave  on  a  Coastline  with  Harbor 

We  shall  follow  a  procedure  of  J.-J.  Lee  (1971)  In  dealing  with  this 
problem.  At  the  mouth  of  the  harbor,  place  a  rigid  fictitious  wall  M  and 
solve  the  problem  in  two  pryts.  First  we  find  the  solution  in  the  region 
outside  the  harbor  with  the  wall  M  in  place,  l.e. , 


SUM 


SuM 


(4.64) 


In  a  simplified  model  of  the  shore  line,  it  might  be  a  straight  line  (see 
Figure  22).  In  this  case,  the  diffracted  wave  would  simply  be  the 


Figure  22  -  Straight-Line  Coast 

plane  wave  reflected  from  the  shore  line.  Note  that  in  this  case,  the 
diffracted  wave  doer  not  satisfy  the  same  radiation  condition  used  earlier. 


an 


However,  with  the  configuration  shown  in  Figure  23,  the  earlier  radiation 
condition  applies.  cannot  really  give  the  solution  outside  the 


Figure  23  -  Circular  Coast 


harbor  for  the  effect  of  the  harbor  has  been  neglected.  The  incident  wave 
<j>^  will  certainly  excite  some  sort  of  motion  inside  the  harbor  and,  in 
turn,  this  motion  will  excite  some  motion  outside  the  harbor  that  will 
radiate  away  from  the  mouth.  Hence  the  motion  in  the  exterior  will  be 
represented  by 


ext 


*I  + 


R 


where  <j>D  satisfies  the  radiation  condition  (4.58).  For  a  point  P  in  the 

K 

exterior,  we  then  know  that 

*R  -  f  J  ‘*RVH01>  '  *R»0v1)]dS 

M  US 

If  we  let  P  be  on  the  boundary  and  remember  that  <t>„  I  _ 

Kn  o 

this  as  follows: 

+R<P>  +  7  j  Vo‘1>dS  -  7  J  *Rv"o1>dS  M-66> 

MUS  M 

We  shall  denote  the  velocity  potential  inside  the  harbor  by  <t>^.  By 
the  Green  theorem,  we  have  for  any  point  P  inside  the  harbor 


(4.65) 

■  0,  we  may  write 
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(1)  A  „  CD 


Vp)  '  J  J  ‘Vo'*'  '  Voi*' lds 


(4.67) 


HuH 


where  we  heve  used  the  Green  function  thet  eetlsfi.ee  the  radiation  con¬ 
dition  because  we  want  to  natch  with  Since  4^1^  “  0,  pert  of 
the  Integral  vanishes.  Letting  P  be  on  the  boundary ,  we  find 

(4.68) 


vp> + 1 ! 


-  i  J 


rH  0V 


MuH  M 

Our  matching  conditions  are  the  following: 


d>  .  m  <b-  +  »  6„  on  M 

Yext  TI  TD  R  rH 


<b  _  ■  on  M 

Text  n  Rn  THn 


(4,69) 


(since  the  positive  direction  of  n  is  opposite  for  <t>R  and  4>H) 

Before  applying  these  in  the  above  equation,  it  will  be  convenient  to 
introduce  some  new  notation.  Let  us  define  the  following  operations: 


I1>  -  * 

H0(1)\  »  \  J  *  Hq(1)  dS 

A 

C>A*  •  i  J  *  Vi1’ dS 


(4.70) 


Then  :he  integral  equations  (4.66)  and  (4.58)  can  be  written  as  follows: 


(I  +  H  (D**uS.  .  ■  H 

U  nQv  ;q>R  nQ  ^ 

(I  +  H  <1>MuHU  -  H  (1)M<b 
KL  *  n0v  'VH  n0  9Hn 


(4.71) 


82 


rmwryjwniBfmi 


Let  ua  denote  the  values  of  <4^  end  <4^  on  M  by 


♦to(P)  -  -  ^(P)  =  F(P),  PCM 


(4.72) 


Then  for  a  point  P  of  M,  we  may  invert  the  above  equations  and  write 


♦r--  <i+^1),,US)'V1)Mf 


♦  -  (I  +  H  U>mu«  -ih  (1)M 
Tl  '  Ov  ‘  0 


(♦.73) 


Since  <4^  -  (jij  +  4>D  +  $R  on  M,  we  have 


(I  +  Ho<1)MuH)-\<1)h  f  -  -  <1  +  C)MuSr\(i)M  »  +  ♦,♦♦* 


or 


-1 


Hq(1)M  F  -  [(I  +  +  (I  +  Hq<1)MuS)'1]  (<))i  +  4,d)  (4.74) 

(DM 


This  gives  F  in  terms  of  <4T  +  <4n.  We  could  proceed  to  find  F,  but 

u  i  U  QNu 

this  is  not  necessary,  for  having  found  Hv  '  F,  we  can  now  use 
Equations  (4.70)  to  find  <4R  or  4>H  on  MuS  or  MuH,  respectively.  Once  we 
know  <j>R  and  <t>R  on  these  boundaries,  we  can  find  them  at  any  point  of  their 


respective  domains.  For  example,  to  find  <4  (P)  for  any  P  inside  the 

H 


harbor,  we  use 


(l)MUH.  .  1  »  (1)M  _ 

VP)  "  -  J  HQv  *H  +  2  tf0  F 


(4.75) 


The  operator  notatio.  '.as  been  introduced  for  the  following  reason. 

If  ve  are  treating  the  exact  problem,  then  these  are  integral  operators 
acting  on  continuous  functions.  However,  if  we  discretize  the  problem  by 
looking  for  say,  only  at  discrete  points  P^  ?2,  .  .  .on  the  boundary. 
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we  then  replace  the  integral  operator  by  a  Matrix  operator  and  the  con- 
tinuoue  function  by  a  vector.  Nevertheless,  the  equationa  in  operator 
fora  remain  the  turn.  Thus  uae  of  the  operator  notation  introduces  a 
certain  economy  in  thought. 

The  net hod  has  been  successfully  tried  by  Lee  (1971)  and  in  a  aore 
elaborate  version  by  Lee  and  Rnichlen  (1972). 

AN  INITIAL-VALUE  PROBLEM  FOR  A  WAVEMAKER 

The  last  example  that  we  shall  consider  differs  from  the  preceding 
ones  in  two  respects:  we  use  a  time-dependent  Green  function  and  we  are 
able  to  give  an  explicit  solution.  The  problem  under  consideration  is 
nothing  but  the  initial-value  problem  for  a  wavemaker  in  a  wall.  It  is 
thus  closely  related  to  the  problem  considered  on  page  26. 

We  suppose  that  some  section  of  the  wall  x  ■  0  is  flexible  and  moves 
in  some  predetermined  way 


x  -  F(yi*»t)  ,  (y,*)eSQ 

We  may  take  the  bottom  to  be  given  by 

y  -  -  H(x,t) 


(4.76) 


(4.77) 


but  eventually  we  take  the  fluid  infinitely  deep  in  order  to  give  an  ex¬ 
plicit  solution.  The  appropriate  boundary  and  initial  conditions  for  the 
problem  are: 


♦x(0,y,*,t)  - 


Ft(y»*»t)  ,  (y,*)eSQ 


0  ,  (y,t)tSQ 


♦n(*.  -  H(x,*),*,t)  -  0  ,  h  •  H(x,x) 


(4.78) 


11a  *  -  0  , 

y 


h  •  «° 


(4.78) 

(Cont'd) 


*tt(x,0,s,t)  +  g*y  -  0 
♦(x,y,*,0)  -  *t(x, 0,*,0)  -  0 
I ♦ | »  | ♦t I »  | | ,  and  | |  are  bounded 

Since  Equation  (4.6)  did  not  involve  time,  it  may  still  be  applied  to 
this  problem.  However,  we  wish  to  make  an  efficient  choice  of  the  Green 
function.  We  shall  suppose  that  it  is  possible  to  construct  a  function 
G  ■  1/r  +  H(P,Q,t)  satisfying  the  following  conditions: 

G^(P;0,n,^;t)  -  0 

'Gv(P;£,  -  H(5,5),C;t)  -  0  if  h  -  H(x,z) 

lim  G  -  0,  if  h  -  <» 

In-  n 

0ttCP!C,0.CiO  +  8cn  -  0  (*.w> 

G(P;Qit)  -  G(P;Q;  -  t) 

« 

or  G  (P;Q}0)  ■  0 

t 

G(P;S»0,C»0)  «  0 

G  -  0(R~2) ,  Gr  -  0(R'3),  Gy  -  0(R-3) 

We  shall  now  apply  (4.6)  to  the  fluid  bounded  by  a  large  vertical 
circular  cylinder  ER  of  radius  r,  and  the  parts  of  the  wall  W,  the  free 
surface  F,  and  the  bottom  B  included  inside  ZR.  However,  we  shall  apply 
(4.6)  to  rather  than  ♦.  Furthermore,  we  shall  take  the  time  variable  in 
♦  to  be  T  and  in  G  to  be  t  -  T,  This  does  not,  of  course,  Invalidate  the 
use  of  (4.6).  We  then  have  the  following: 

4irt>r(P;T)  ■  f  l-  ♦t^(0,n,c;T)G(P;0,n,i;;t  -  t)  +  ♦^Jdndc 

w 

+  J  l4’tn(t.0;C;T)C(P;C,0,C;t  -  t)  - 


•5 


r.  ii^g^rvr’nvg.r 


where  we  have  already  used  the  boundary  conditions  on  B  and  the  asymptotic 
condition  for  large  R  to  eliminate  the  Integrals  over  B  and  ER.  We  may 
also  exploit  the  boundary  conditions  on  W  and  F  to  write  this  equation 
as  follows: 

4ir$t<P;x)  -  -  J  Ftt(n,C;t)G(P;0,n,i;;t  -  x)dnd£ 


+  |  f  [-  *ttt(e,0,C;T)G(P;C,0,C;t  -  T)  +  *tGttldWc 
JF 


1 


J  Ftt^T)G(t  “  T)dndc 


'  i  /  It  ^tt<T)G(t  -  T)  +  VT)Gt(t  -  'DldWe 

F 


(A. 80) 


J  Ftt'T)G(t  "  T>dndt 


+  I7  J  [*n(T)G(t  -  T)  -  i  *t(T)Gt(t  -  T)]dCdC 


F 

He  now  integrate  this  equation  from  0  to  t: 

t 


4x*(P;t)  -  4ir*(P;0)  -  -  J  dx  J  Ft(;(x)G<t  -  x)dnd£ 


+  J  [*n(t)G(0)  -  i  ♦t(t)Gt(0)]d5dC 


-  J  t*n(0)G(t)  -  i  ♦t(0)Gt(t)JdCdi; 


(A. 81) 
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If  we  now  make  use  of  Che  initial  conditions  for  G  and  of  Equations  (1.9) 
and  (1.10),  we  say  write  this  equation  as  follows: 

4TT<KP;t)  -  4tt*(P;0)  -  J  dT  J  Ftt(t)G(t  -  T)dnd£ 

0  S 

°  (4.82) 

-  J  tYt(C,C,0)G(P;C,0,?;t)  +  Y(S,S,0)Gt(P;5,0,C;t) ]d£dC 
F 


In  fact,  we  have  assumed  that  4> (P ; 0)  “  0  and  ♦  (£,0, C»0)  *  0  as  initial 
conditions  for  <$.  Hence  (4.82)  becomes 


4TT<t»(P,t) 


^t 

J  dT  |  Ftt(n,^,T)G(P;0.n,C;t  -  T)dr)dC 


(4.83) 


This  selves  the  problem  formally,  but  we  are  left  with  the  problem  of 
finding  a  Green  function  satisfying  (4.79).  This  function  can  be  con¬ 
structed  for  either  H(x,z)  ■  h  ■  const  or  h  -  We  give  it  here  only  for 
h  ■  “.  We  give  first  the  result  without  the  condition  G^(P;0,ri,^;t)  «  0: 

00 

G0(P,Q,t)  -  7-  7+  2/  dk  ek(y+n)  JQ(kR)[l  -  cos(/£Tt)] 

0  (4.84) 

00 

-7  +  7  ~  2  J  dk  ek(y+n)  J_(kR)  cos  (v'gkt) 

r  r.  j  u 

0 


where 


r  -  [(x  -  C)2  +  (y  -  n)2  +  (z  -  O2] 

rx  ■  [(x  -  C)2  +  (y  +  n)2  +  (*  -  c)2 
2  2 

R  -  ((x  -  O  +  (*  -  O  ] 


1/2 

1/2 

] 


In  order  to  construct  a  Green  function  that  also  satisfies  G^(P;0,n,4»t)  * 
0,  we  use  the  method  of  images.  Define 
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r2  -  l(x  +  5)  +  (y  -  ri )  +  U  -  5)  ] 


r3  -  [(x  +  O2  +  (y  4  n)2  +  (*  -  C)2] 


1/2 


R2  -  [(x  +  C)2  +  (z  -  C)2] 


1/2 


Then 


G(P;Q;t)  «7  +  -£-  +  —  +  f- 
r  ^  r2  r3 


(4.85; 


w 

-  2  J  dk  ek(y+n)  [JQ(kR)  +  J0(kR2)l  cos  t^J 


If  we  substitute  this  into  (4.83),  we  obtain 


{ 


4ir<P(P;t)  -  -  2  |  dnd?|  - 

'S. 


(y  -  n)2  +  (*  -  c)2 


+  - - - - - ---1  f  f 

S*2  +  (y  +  n)2  +  (z  -  n)2  J -l 


tt 


(n,c;t)dx 


(4.86: 


J 


+  4  I  dndc  J  dk  ek(yfn)J0(ktx2  +  (z  -  O2]  J  F^(n,C;t) 


tt 


so  0 


cos  /gk(t  -  x)dT 


If  we  make  the  special  choice  (2.38), 

F(y,z,t)  -  F(y,z)  sin  at 

the  integrals  with  respect  to  T  can  be  evaluated.  An  interesting  problem 
is  then  to  find  the  asymptotic  form  of  the  solution  as  t  +  *.  This  should 


88 


and  does  agree  with  the  form  obtained  by  the  Fourier-Integral  method.  The 
asymptotic  form  may  be  written  as  follows: 


My+n) 

JQ(kR)  -  2irk0  cos  at  e  Vk0R) 


The  use  of  time-dependent  Green  functions  in  the  manner  shown  above 
is  apparently  first  due  to  Volterra  (1934) .  It  has  been  further  exploited 
by  Finkelstein  (1957),  Wehausen  (1967),  and  W.-C.  Lin  (1966).  Many  people 
prefer  to  work  entirely  with  frequency  by  initially  taking  Fourier  trans¬ 
forms  (or  Laplace  transforms) . 
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PART  5:  THE  METHOD  OF  MULTIPOLES 

The  method  of  multipoles  can  be  used  in  either  two  or  three  dimensions 
and  consists  essentially  in  adding  together  singularities  of  higher  and 
higher  order  with  intensities  to  be  determined  by  the  boundary  conditions. 

In  the  case  discussed  here,  the  singularities  are  chosen  to  satisfy  the 
free-surface  and  radiation  conditions.  The  Intensities  are  then  selected  to 
satisfy  the  boundary  condition  on  an  oscillating  body. 

The  method  occurs  in  various  guises,  but  in  the  type  of  problem  we  con¬ 
sider,  it  was  apparently  first  used  by  Ursell  (1949),  who  treated  the 
heaving  motion  of  a  circular  cylinder.  The  analogous  problem  for  a  heaving 
sphere  was  later  considered  by  Havelock  (1955) .  Others  who  have  used  it  are 
Tasai  (1959),  Porter  (1960),  and  C.M.  Lee  (1968),  but  the  list  is  not 
exhaustive. 

The  present  approach  to  the  method  differs  somewhat  from  the  usual  one, 
but  as  will  be  shown  later,  the  results  are  the  same.  The  present  approach 
seems  to  have  the  advantage  of  showing  from  the  beginning  that  the  sum  of 
singularities  usually  employed  really  is  sufficient  to  solve  the  problem. 
Once  again  the  method  is  illustrated  by  the  problem  of  the  forced  harmonic 
motion  of  a  body  in  a  free  surface.  The  treatment  will  be  two  dimensional 
and  will  have  much  in  common  with  Chat  of  diffraction  from  a  vertical  plate. 

FORCED  HARMONIC  MOTION 

Figure  24  shows  schematically  the  physical  situation.  He  shall 


Figure  24  -  Oscillating  Body  on  the  Free  Surface 
suppose  that  there  is  a  complex  velocity  potential  of  the  form 
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■  *B.wt 


•  •  •  I 


(5.9) 


ft 


g(3)  -  cz  +  aQ 


>  a 


where  c  and  a  are  real.  We  may  also  conclude  from  (5.4),  just  as  in  the 
n 

diffraction  problem,  that  c  *  0.  Furthermore,  a^  may  be  absorbed  into  the 
stream  function  without  loss  of  generality,  so  that  we  finally  have 


CO 

g(a)  -  X  -J  ,  |*|  >  a  (5.10) 

n-1  zn 

f 

\  One  of  the  important  properties  of  power  series  representing  analytic 

i  functions  is  that  they  converge  absolutely  up  to  the  nearest  singularity. 

Since  the  function  g  must  be  analytic  everywhere  outside  the  body,  we  know 
{•  that  (5.10)  must  converge  at  least  up  to  a  circle  just  containing  the  body, 

I  as  shown  in  an  example  given  in  Figure  25.  If  the  body  is  a  semicircle, 


i 


l 


SERIES  CONVERGENT 
OUTRIDE  THIS  CIRCLE 


Figure  25  ~  Region  of  Convergence 

then  the  series  (5.10)  will  converge  right  up  to  the  body.  Moreover,  since 
there  is  no  singular  behavior  of  g  on  the  circle  representing  the  body,  the 
radius  of  convergence  must  be  smaller  than  the  radius  of  the  body  itself. 

It  will  also  be  true  that  the  radius  of  convergence  of  (5.10)  for  a  non¬ 
circular  body  with  smooth  boundary  will  be  smaller  than  the  circle  just  em¬ 
bracing  the  body.  However,  this  is  no  guarantee  that  it  will  hold  right  up 
to  the  boundary  of  the  body.  Thus,  at  the  moment,  we  seem  to  be  con¬ 
strained  to  treat  only  a  body  of  circular  section  with  center  in  the  free 
surface.  For  the  time  being  we  shall  do  this  and  later  consider  what  to 

do  if  the  cylindrical  section  is  not  circular. 

i 

i 

I 

I 


Now  each  coefficient  after  the  first  is  purely  imaginary,  and  it  is  possible 
to  show  that  each  term  individually  satifies  the  tree-surface  condition 
and  radiation  condition.  Define 


«„(«> 


_j_  +  i  i- 

n+1  n  n 

2  Z 


Then  it  is  easy  to  confirm  that 


(5.18) 


* 

Im{m  (xHO)  v  ikAm  }  ■  0 
n  On 


Hence,  b  m  (z),  b  real,  also  satisfies  this, 
n  n  n 

We  assert  that  we  may  write 


Va  f 
n  n 


Vi  + 


b  m  (z) , 
n  n  '* 


n*2 


»^2 


3* 


real  (5.19) 


of  course,  the  manipulations  leading  to  (5.17)  do  not  prove  this,  but  once 
the  combination  defining  has  been  recognized ,  it  is  not  difficult  to 
prove  (5.19).  For  example,  we  may  do  this  by  using  the  relationship 

fn  '  -  72  ["n  +  <n+1>  £n«] 

*0 


Let  us  now  consider 
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*!-• 


-lk„i  '  lk0C 


I  hr*1 


-J 


*  ik0(C-*» 


d( 


W*  shall  make  the  change  of  variables 

kQ((-*)  -  -ks  or  k  -  kQ  ^  ,  ( 


(5.20) 


Formally,  this  gives  Immediately 


•J 


0  -Iks 
e _ 

k-k-. 


dk 


dk 


where  the.  path  of  Integration  Is  still  to  be  determined.  Let  us  take  our 
path  In  the  (-plane  as  shown  In  Figure  27  (two  cases  are  shown). 


Figure  2 


Paths  of  Integration  in  (-Plane 


The  transformed  paths  in  tha  k-plane  will  than  appaar  as  shown  in  Figure  28 
Tha  naw  paths  ara  a  rotation  and  otratching  of  tha  original  paths. 


k-PLANF. 


Figure  28  -  Transformed  Paths  in  the  It- Plane 


Let  us  now  examine  the  arcs  joining  tha  path  to  «°  with  the 


positive  x-axis  (see  F ignis  29) . 


Figure  29  -  Arcs  £_  Joining  the  Path  to 


W«  find  easily  that  ths  contribution  of  goes  to  tsro  as  ?.  *  “•  Hanca 
wa  aay  deform  the  path  to  ona  slong  tha  x-axis  except  for  dodging  around 
kQ  (sea  Figure  30). 


Figure  30  -  Path  along  the  x-Axis 


If  we  choose  the  path  around  kg  as  a  semicircle  of  radius  e  and  then  take 
the  limit  as  e  -*■  0,  we  will  obtain  a  Cauchy  principal- value  integral  plus 
a  half  residue: 


f^*) 


f  «rikl  “ik08 

«  7—r —  dk  -  irie 

J  k’ko 


(5.21) 


Now  consider  the  radiation  condition  (5.3).  In  order  to  satisfy  this, 
we  must  find  the  asymptotic  behavior  of  f  as  x  +  i  J0.  If  we  use  the  series 
(5.19)  as  part  of  the  representation,  it  is  obvious  that  we  may  discard  the 
contribution  of  the  terms  in  n»n  as  x  •+■  ±  °°.  Hence  we  need  to  examine  the 
asymptotic  behavior  only  of 

-ik-s 

a^s)  +  bce  °  (5.22) 

It  is  evident  that 

~ikn* 

f(e)  — •  bge  as  x  -v  +  00  (5.23) 

If  we  write  (5.22)  in  the  form 
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As  in  the  case  discussed  on  page  50,  where  the  path  in  Figure  11  has  been 
replaced  by  that  in  Figure  12,  we  may  replace  the  integral  from  to  f» 
by  the  closed  path  shown  in  Figure  31. 


Figure  31  -  Equivalent  Path  of  Integration  for  f(s) 
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We  now  apply  the  radiation  condition  (5.3)  to  f  and  f^,  adding  appropriate 
subscripts  to  the  coefficients: 


fc'  +  k0£s'  ~  -lk0(bc0  +  lbc0  >  +  k0(bs0  +  lbs0  >  ■  0  aa  x  *  +  “ 

fc’  -  k0£.’  ~ -V^Vl  +  bc0  +  ibco'> 

-ko(-2,rlSi  +  bsi  +  ib9o')  *°  aax 


These  equations  yield 


t  •« 

^sO  ""  ^cO 


■«d»  b 


cO  "  bs0  ‘“si 


Let  us  now  go  back  to  (5.21)  and  substitute  this  form  for  f^: 


f  *  a 
c  cl 


l  k'k« 


dk  ~  trie 


-ikQz 


+  '''“.I  +  iad)a 


-lk02 


+  7  ■  D  IE  | 

n«2  cn  n 


£„  "  aal 


;f 


-ikz 


k-k„ 


dk  -  Trie 


-ikQz 


+  Tr(-aci  +  iasl)e 


-ik0z 


+  T  b  m 
sn  n 


n*2 


(5.25) 


These  simplify  to  the  following  forms : 


‘cl  J  k-k, 


«  “ikOz  v- 

-—  dk  +  rca  .  e  +  X,  b  m 

,  si  *-«  cn  n 


(5.26) 


^s  asl  J  k-k. 


'lk02  v~ 

dk  -  iracl  e  +  ^  banm0 


Before  returning  to  the  problem  of  satisfying  the  boundary  condition 
on  the  body,  we  make  one  more  digression,  namely,  to  explain  why  this  is 
known  as  che  method  of  multipoles.  Suppose  that  we  wish  to  construct  a 
function  behaving  like  a  source  of  strength  Q  cos  at  and  vortex  of  intensity 
T  cos  at  at  the  point  c  in  the  lower  half-plane,  i.e.,  like 


T+iQ 


log  (z-c)  cos  at 


and  also  satisfying  the  free-surfaoe  and  radiation  conditions  and 
vanishing  as  y  +  -  «.  This  problem  can  be  solved  (e.g.,  by  the  reduction 
method)  and  we  find 

G(z,c)  -  log  (z-c)  +  log  (z-c) 

- 1  ^  f  e  k-v"  - dk] cos  c,t 


(5.27) 


-ik-(z-c) 

-  i(r-iQ)e  sin  at 


■  G  cos  at  +  G  sin  at 
c  s 


The  method  of  multipoles,  as  applied  to  the  problem  we  have  been  con¬ 
sidering,  would  consist  in  assuming  that  the  solution  can  be  expressed  as 
a  sum  of  G  and  its  derivatives  evaluated  at  c  ■  0,  i.e.,  that 


f  -  anG  (z,0)  +  a..G  (z,0)  +  a«G  (ffjO)  +  ...,  a  real 

C  U  C  X  C  Z  C  Q 

f  “  a»G  (z,0) 


(5.28) 


Since  there  is  no  circulation  in  the  problem  we  are  considering,  we  may  set 
r-0.  Furthermore,  there  is  no  loss  in  generality  in  setting  Q*1  since  each 
summand  is  multiplied  by  a  real  coefficient  to  be  determined  later.  If  we 
also  set  c“0,  (S.27)  then  takes  the  following  form: 


Gc(z,0) 

Gs(z,0) 


dk 


\5.29) 


Now  if  we  examine  the  terms  of  (5.26)  before  the  summation,  we  see  that 

by  setting  agj*0,  bsn*0,  we  have  exactly  the  same  series  as  (5.28). 

Although  (5.26)  may  appear  to  be  more  general  because  of  the  presence  of 

the  terms  with  a  .,  and  b  ,  this  is  only  an  apparent  generality,  an 

appropriate  shift  of  the  time  axis  can  achieve  a8i“^8n“®*  Hence  the  method 

we  have  developed  is  equivalent  to  the  method  of  multipoles. 

Next  we  must  satisfy  the  boundary  condition  (5.7)  by  proper  choice  of 

the  coefficients  a  .,  a  .,  b  ,  and  b  .  For  the  case  to  which  we  are  pre- 
ox  si  cn  8ii 

sently  restricting  ourselves,  that  of  the  semisubmerged  circular  cylinder, 
this  condition  may  be  given  a  slightly  simpler  form.  In  order  to  exploit 
symmetry,  we  Introduce  the  angle  y  shown  in  Figure  32. 


Figure  32  -  Definition  of  Angle  y 


As  in  (4.13),  we  shall  write 


Vn(x,y,t) 


V  coa  at  +  V  sin  at 
nc  ns 


(5,30) 


Then  (5.7)  takes  the  form 


A  m  -A  sin  Y  i  4  cos  Y  «  V  (y) 
Ycn  Ycx  '  Ycv  1  nc  1 


b  -  -A  sin  y  +  v  cos  y  *  v  (Y) 
rsn  Ysx  1  sy  1  ns 


(5.31) 


It  will  be  sufficient  to  discuss  one  of  these  although,  as  we  shall  see, 
the  two  must  be  solved  together.  From  the  solution  (5.26)  we  may  compute 


f  **  +  ia  , 
c  cl 


w 

f  fc 


■: - kdk  -  lira  ,k„e 

kQ  si  0 


+  ^  b  m 


(5.32) 


From  (5.18) 


M  •  -(n+1>  -KS?  lk0  ~S+T  "  (5-33) 

z  z 

iv 

Into  (5.32)  and  (5.33),  we  substitute  z  ■  -ire  and  separate  real  and 
imaginary  parts.  This  somewhat  tedious  computation  yields  (perhaps)  the 
following: 

OO 

r  -krcos  y  .  , 

*  -  .  ,  t  2 - Wtkr.ln.Jl) 

Ycx  cl  J  k-kn 


-kgrcos  y 

-  ^g^k^e  8in(k^rsin  y) 


OO  — 

+  S  bc,2p  ^2p+2  ("1)P  LC0S(2P+2)Y  "  2?+r  8i*<2P+1>Yj 


(5.34) 


+  Z  bc 1 2q+i  (“1>q+1  [8in(2qf3)Y  +  2i?2  C08<2q+2>Y] 
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<p 

VJ 


"*aci 


i 


-krcos  Y  ,,  .  x 

e _ ’costkrsin  y) 


k-k. 


kdk 


-kgrcosy 

+  ira^e  cos(kgrsinY) 


V 


+  ^  bc,2p  ^2  ("1)P+1  [sin<2P+2)Y  +  "2p+i  cos(2p+l)Y 


] 


+  ^  bc,2q+l  r2q+3 


[• 


1c  r 

(-l)q+1  |co8(2q+3)Y  -  o?Zo  sin(2q+2)Y 


2q+2 


It  seems  reasonable  now  to  skip  over  some  details  and  to  observe  that 
(5.31)  takes  the  following  form: 


aclVY>  +  °siMcW  +  ^  »cn»cCv)  -  vncCr> 


n=»2 


“clLs(Y>  +  SlVY>  +  ^  »en»8(V)  ■  Vm(Y) 

n»2 


(5.35) 


There  are  evidently  an  Infinite  number  of  unknown  coefficients,  but  also, 
of  course,  an  infinite  number  of  values  of  y  where  these  equations  must  be 
satisfied.  We  shall  not  discuss  the  methods  that  have  been  used  to  solve 
approximately  these  equations  by  truncation.  Unfortunately,  the  various 
functions  have  no  orthogonality  properties. 

We  must  now  consider  the  question  of  what  to  do  ii  the  cross  section 
of  the  body  is  not  a  semicircle,  for  as  we  have  seen  earlier  (see  page  97), 
we  have  no  reason  to  believe  that  the  circle  of  convergence  of  the  function 
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g(z)  in  (5.10)  will  be  inside  the  body  boundary  (plus  its  reflection). 
What  we  do  in  this  case  is  to  make  a  preliminary  mapping  of  the  exterior 
of  the  body  plus  its  reflection  onto  the  exterior  of  a  circle.  However, 
the  mapping  chosen  imposes  a  restraint  upon  the  body:  it  must  Intersect 
the  free  surface  perpendicularly. 

Let  C  be  the  plane  of  the  circle  and  z  the  physical  plane  containing 
the  body  in  question  (see  Figure  33) . 


Figure  33  -  t,-  and  z-Planes 

Then  it  is  known  that  we  can  map  the  exterior  of  the  circle  onto  the 
exterior  of  the  body  and  the  £-axis  into  the  x-axis  by  a  mapping  of  the  form 


00 

_  .  r 

z  -  5  +  ^  ~  , 

C  real 
m 

(5.36) 

m*l  C 

Let  us  suppose  that  f(z)  is  one  of  f^  or 

f  .  Then 
s 

f(z)  -  f(z(C)) 

=  F(5) 

(5.37) 

defines  an  analytic  function  in  the  £-plane  that  can  be  considered  as  a 
velocity  potential  there.  What  happens  to  the  free-surface  condition? 
Since  F* (C)  ■  f * (z)  z'(C)»  we  find 

I 

f'(z)  +  iknf(z)  -  +  ik  f(C) 

0  z  (?)  0 

Hence  condition  (5.2)  becomes 
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or,  since 


is  real. 


I  jFt(i+iQ)  +  ik  F(C+iO)l  *  0 
mU(C+iO)  °  ) 


z  (5+iO)  *  1  -  Z  mC  £ 


00 

Im  |f'(£H-0)  +  ik0  jl  -  Z  mCm5'lr'1j  F(fc,+iO)|  -  0  (5.38) 

This  is  the  boundary  condition  that  must  be  satisfied  on  the  real  axis  out¬ 
side  the  circle. 

We  shall  now  proceed  similarly  to  the  eSrlier  case.  Define 


00 

Gtt)  -  f'(C)  +  ik0  [l  -  Z  ~JlJ  F(C) 


(5.39) 


This  function  can  be  extended  by  reflection  into  the  whole  complex  plane 
and  furthermore  can  be  shown  by  easy  arguments  to  be  bounded  outside  some 
circle  |c|  *  a  containing  the  mapping  circle.  As  before,  we  find 


-5. 


n*l  C 


-  ,  an  real,  U|  >  a 


(5.40) 


But  then,  since  G(£)  must  be  analytic  right  up  to  the  boundary  of  the 
mapping  circle,  the  Laurent  series  (5.40)  must  converge  right  up  to  and 
including  this  circle. 

We  are  now  left  with  the  problem  of  integrating  the  differential 
equation 


00 

i  r  .  _ ,  mC  t  m  »  a 

'  tt> +  lko  [  -I^r] -  2  7 


(5.41) 
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The  integral  can  be  found  In  textbooks  on  differential  equations  (e.g., 
Kanke,  "Different ialgleichungen:  Losungsmethoden  und  Losungen,"  p.  16) 
and  may  be  written  as  follows: 

F(C)  -  exp{-ik0U+£CnC'm)}  |b  +  £  «n  JC  C“a  exp{ik0tc+rxac'‘m]}6c| 

'  11  *  (5.42) 

where  the  path  of  integration  is  below  the  mapping  circle  and  B  is  complex. 
Let  us  define 

Fn(C)  -  expt-il^tC+IC^'®)}  c“nexp{ik0[C+ZCmrm)}dC  (5.43) 


Then  (5.42)  may  be  written  in  the  form 


CO 

F(C)  -  B  exp{-ik()[C+£CinC“m]>  +  ^  »n  Fn(C) 


(5.44) 


If  we  integrate  F  once  by  parts  for  n  >  2,  we  find 
n  ~ 


-n+1 

p„  ■  W  +  •pp{-lko«'+!:c.'  »  11 


f  **0  [l  "X  m+ll  D-l  eXP^lk0li:+!:CBC  lJdl: 

00  s  J  s 

■  r^T  n-1  +  n-1  ik0Fn-l  "  n^l  mCmFm+n 


Let  us  replace  n  by  n+1  and  rearrange  the  result  as  follows: 


Fn(?)  “  ikQ  Fn+1  ”  X  mCmFm+n+l  “  ikQ  ^n 

.n-1 


(5.45) 


n  *  1)2|3|  .... 
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If  possible,  ve  should  like  to  do  something  like  we  did  on  pages 
99-100,  i.e. ,  to  replace  the  awkward  integrals  ...by  polynomials, 

something  analogous  to  the  mn(s)  of  (5.. 18).  Equation  (5. AS)  is  the  analog 
of  (S.13).  However,  to  manipulate  terms  as  we  did  in  (S.16)  and  (5.17) 
would  seem  to  require  a  remarkable  insight.  Let  us  write  out  (5.A5)  for 
n  »  2,3,  ....  for  several  terms: 


F2  -  'jjj-  F3  -  -  2C2  P5  -  3C3  Fg 


1  1 
lko  ;2 


F3  -  F4  -  C1  F5  -  2C2  F6 . 


Fj  “  ,  Fg  “  C-  Fj.  "*eee< 

A  ik„  S  16 


1  1 
lkoc3 


1  1 
ik0C4 


(5.46) 


F5  “  ik^  F6  “*•* 


1  1 
ik0  C5 


It  seems  evident  that  Fft  must  be  expressible  as  u  series  starting  with  C~u: 


*  s  <n> 

r  -ilV 

n  ?n  ££  ^k 


(5.47) 


If  we  substitute  this  expression  into  the  defining  equation  for  F 


*n'+ivjl- .fSjrll  -i- 
”  0  L  ±1  "  c“ 


it  is  not  difficult  to  establish  the  following: 


n  •  -T2L  B  -  Hti.  b  +  r  b 

Br  ikQ  *  B1  k  2  »  B2  ikQ  B1  +  C1B0 


and  in  general  that 
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«-*>  vA 


(3.48) 


(n) 


> 


<1>1 


Although  n  could  proceed  with  this  exprescion  for  F^,  It  would  require 
cereful  attention  to  detail  in  separating  real  and  Imaginary  parte  In  the 
coefficients  B  q>2. 

In  order  to  avoid  this  difficulty,  it  has  been  cuetomary  to  proceed 

somewhat  differently.  Instead  of  the  sum  £  a  F  in  (5.44),  we  shall  try 

n  n 

to  replace  it  by  another  sum  £  d  M  ,d  real,  where  the  M  will  be  defined 
7  nun  n 

below.  We  shall  again  suppose  that  has  the  form 


M 


n 


i-I 

Cn  k-0 


(5.49) 


but  we  shell  begin  by  imposing  only  the  condition 


Im 


[Mn  +  ik0 


when  H  -  0 


(5.50) 


Of  course,  F  satisfies  this  condition  also,  so  that  we  must  expect  to 
n 

impose  some  additional  conditions.  Substituting  (5.49)  into  (5.50),  we 
obtain  the  following  equation  (we  temporarily  drop  the  superscript  u) : 

{‘Vc  +  «V!  '  "V  '"l  + 

00 

+  2^  £ik0bk+l  -  (k+n)bk 
*  • « 

If  we  now  set  bk*bk  +ib^  and  C*f;+iO,  the  coefficients  of  each  power  of 
£  yield  the  following  equations: 


k-1 

-  1*0  Z  <k-r> c. 


r**0 


-k-1 


l 


(5.51) 
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W«  have  not  yet  shown  that  5]  a  F_  can  b.-»  teplaced  by  a,F.+  ^  d  H 

n-1  nn  1  1  n-2  n  R 

with  dQ  real.  With  son*  further  wanipuiation,  we  can  establish  the  follow¬ 
ing  relation  between  F  and  H  : 

n  n 

\  -  -  TT  +  k  (k  *  Sa  )  — *]  V«  +  V  (k  +  S3  )  Fn« 

+  k02  \  (<r-l>  (S  +  Vr-k-ij  (5.55) 

r“4  l  k-1 

Since  the  coefficients  on  the  right-hand  side  are  all  real,  it  is  evident 


that  a  sum  Z*  d  M  with  d  real  is  equivalent  to  a  sum  A  a  F  with  a. 

„  n  n  n  ,  n  n  n 

n-2  n-2 

real.  Since  the  coefficients  have  not  yet  been  determined,  it  is  not 

necessary  to  know  the  relation  between  the  {a  }  and  the  {d  }.  Instead,  we 

n  n 

shall  determine  the  d  directly. 

n 

Equation  (5.44)  for  F(£)  now  takes  the  form 

00 

F(C)  -  B  exp{-ikn[t+EC  ?“”*]}  +  a.F.  +  X  d  1 M  ,  d  real 

0  m  1  1  n-2  n  n  n  (5.56) 

The  radiation  condition  and  the  boundary  condition  on  the  body  remain  to  be 
considered.  For  this  purpose  there  may  be  some  gain  in  replacing  (O  by 
f^(*(C))  where  f^  is  defined  by  (5.12)  with  n-1.  This  is  not,  of  course, 
the  same  function  as  F^(C).  but  it  can  be  used  in  (5.56),  as  we  shall  see 
below: 


a  sum  ^  a 
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C  i  V^2- 

f  «xp  ik0tc+£Cac“*l 

j  r  C  ^ 


•xp  ik0u+£cttrBi 


'-U*±x\« 

L  n“l  C  J 


where  eR  in  reel.  Evidently  the  difference  between  f^(s(C))  end  F^(£)  will 

QO 

y 

be  a  sua  enF  which  in  turn  can  be  replaced  by  a  sum  of  M  . 

n«2  n 

Hence  we  may  replace  (5. S3)  by 

F(C)  -  B  exp{-lk0tC+ECmC“®]}  +  a^UU))  + 


+  Z-*  d  M  ,  d  real 
,  n  n’  n 
n-2 


(5.57) 


where  these  d  are  different  In  general  from  those  of  (5.56). 
n 

With  this  change ,  we  may  take  over  completely  our  earlier  relations 
on  page  104  that  were  derived  to  satisfy  the  radiation  condition,  for  the 
Mq  play  no  role  In  this  condition. 

In  order  to  satisfy  the  condition  (5.7)  on  the  body  in  the  physical 
plane,  we  must  establish  a  relation  between  this  condition  and  one  in  the 
C-pltme.  Let  us  first  rewrite  (5.7)  for  the  case  where  the  parameter  is 
no  longer  arc  length.  Let  it  be  p.  then  (5.7)  becomes 


Im  f  (s(p) )  7-7 


-xM - - 


(5.58) 


I*  (P)| 


Let  p  now  be  a  parameter  in  the  C~ plane  describing  the  mapping  circle: 
C(p).  Then  the  body  in  the  physical  plane  is  described  by 
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s(p)  -  «(C(p>) 


and 

Sine* 

we  find 


«  (p)  •  «  (C(p))  C  (p) 

Ms)  «'(0  -  f’(C) 

t  (s(p))  t'(p)  -  f'(  *)  *'«>  C*  <P)  -  MO  Mp> 


condition  (5<S9)  then  becomes  the  following  itv  the  IJ— planet 


(5.59) 


To  develop  the  analogs  of  (5.35)  would  involve  us  in  more  details 
than  seem  appropriate  to  carry  through  in  these  lectures.  We  terminate 
the  discussion  of  this  method  with  the  remark  that  it  is  possible  to 
develop  this  into  a  computational  technique  for  bodies  intersecting  the 
free  surface  perpendicularly,  as  has  been  rhown  by  Tasai  (1959)  and  Porter 
(1960) . 
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